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1. Introduction. This is a review paper of the Bukhgeim-Klibanov method (BK). The author considers 
this paper as an introductory material for BK. Because of many publications on BK, the author restricts 
himself to citations of those works which he is best familiar with. An interested reader can find further 
citations as well as analytical details in publications cited here. Three topics are discussed in this paper: (1) 
Global uniqueness for Multidimensional Coefficient Inverse Problems (MCIPs) with single measurement data 
via BK, (2) Global stability both for MCIPs and some ill-posed Cauchy problems via Carleman estimates, 
and (3) Related convergent numerical methods both for MCIPs and some ill-posed Cauchy problems. 

BK was introduced in three originating papers of Bukhgeim and Klibanov in 1981 [3S1 [73]. Until 
now BK remains the only technique enabling to prove global uniqueness and stability theorems for MCIPs 
with single measurement data. The state of the art in the field of Inverse Problems in 1981 is well reflected 
in the following citation from the paper j35j . '^Uniqueness theorems for multidimensional inverse problems 
have at present been obtained mainly in classes of piecewise analytic functions and similar classes or lo- 
cally... Moreover, the technique of investigating these problems has, as a rule depended in an essential way on 
the type of the differential equation. In this note a new method of investigating inverse problems is proposed 
that is based on weighted a priori estimates. This method makes it possible to consider in a unified way a 
broad class of inverse problems for those equations Pu = / for which the solution of the Cauchy problem 
admits a Carleman estimate... The theorems of %1 were proved by M.V. Klibanov and those of ^2 by A.L. 
Buhgeim. They were obtained simultaneously and independently.^^ 

While the paper [35] contains announcement of results and an indication of the proof, papers [36] [74] 
contain first complete proofs. Also, see, e.g. [37l[48l[ni[l5l[76l[7a[78l[8a[79l[84l[8a Sections 
1.10 and 1.11 of the book [14] for some follow up publications of these authors on BK. Prior publications 
[35] [36] [74] only the so-called local uniqueness theorems were known for MCIPs with single measurement 
data. The term local means here that unknown spatially dependent coefficients were assumed to be either 
piecewise analytic functions, or functions represented via truncated Fouricr-likc series, or sufficiently small 
perturbations of constants. At that time the common desire of many mathematicians working on MCIPs 
was to prove global uniqueness theorems. That is, to prove uniqueness for the case when the unknown 
coefficient a{x) ,x G M",n > 2 satisfies only some natural conditions, such as, e.g. a € C'^. However, it was 
unclear at that time how to do this. Indeed, standard methods were basically based on integral equations 
and did not work for this goal. The absence of global uniqueness theorems was the main stumbling block 
in the field of Inverse Problems in 1970s. This is what has originally motivated the author in 1979 to think 
about moving away from traditional techniques. The author has spent two years to figure out the solution. 

This paper is focused only on MCIPs with single measurement data. BK is based on a special use of 
Carleman estimates for MCIPs. Roughly speaking, as soon as Carleman estimate is valid for a PDE operator, 
BK can be applied. Therefore, the generality of BK is due to the fact that Carleman estimates are valid for 
hyperbolic, parabolic, elliptic, the non-stationary Schrodinger and some other operators. Schematically, BK 
consists of the following two steps: 

Step 1. Given an MCIP, figure out whether a proper Carleman estimate is valid for the PDE operator 
of this problem. If not, derive a proper Carleman estimate (if possible). 

Step 2. Given a proper Carleman estimate, apply BK. 

Carleman estimates were originated in 1939 in a remarkable work of a distinguished Swedish mathe- 
matician Torsten Carleman [42j . Since then and up to now they have been traditionally applied by many 
authors to proofs of uniqueness theorems for various ill-posed problems for PDEs with the Cauchy data 
on non-characteristic hypersurfaces, see, e.g. the paper of Calderon [40] and the book of Hormander [57] . 
While in these references Carleman estimates were used for functions with compact support, the book of 
Lavrcnt'ev, Romanov and Shishatskii [112] (Chapter 4) uses functions with non-compact support. As a result, 
the technique of this book allowed to prove not only uniqueness but Holder stability results as well for some 
ill-posed Cauchy problems. Thus, it was briefly noticed in earlier publications jTS] [HI] with the reference to 
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[112] that the applicability of BK to an MCIP usually automatically implies the Holder stability estimate 
for this MCIP. Since the author is concerned with MCIPs, multiple quite interesting works on applications 
of Carleman estimates to various ill-posed Cauchy problems are not cited here. 

BK applies Carleman estimates, in a specially designed way, to proofs of global uniqueness and stability 
results for MCIPs. MCIPs are substantially different from those Cauchy problems. Indeed, in such a Cauchy 
problem all coefficients of the corresponding PDE operator are known, the initial condition is unknown, 
and the Cauchy data at a non-characteristic hypersurfacc are known. It is required then to reconstruct the 
solution u of the corresponding PDE. On the other hand, in an MCIP at least one of coefficients, a {x) , of 
the corresponding PDE operator is unknown, the initial condition is known, and the lateral Cauchy data 
are known as well. It is required then to reconstruct the pair of functions (a, u) . Those Cauchy problem are 
linear and any MCIP is nonlinear. 

The term "MCIP with single measurement data" means the problem of the recovery of one of coefficients 
of a PDE from a boundary measurement generated by a single set of initial conditions. In the case of either 
the point source or a plane wave this means either a single position of that source or a single direction of 
that incident plane wave. More generally, this is a single pair of initial conditions for a hyperbolic PDE 
and a single initial condition for a parabolic PDE. Sometimes a few initial sets of initial conditions are 
allowed. In the case when k coefficients are unknown, k sets of initial conditions are allowed, which means k 
measurements. MCIPs with single measurement are non-overdetermined ones. The non-overdetermination 
means that the number of free variables in the data equals the number of free variables in the unknown 
coefficient. The single measurement case is the one with the minimal amount of the available information. 
Hence, this is the most economical way of data collection. In particular, in military applications the single 
measurement case is far preferable to the case of many measurements. This is because an installation of 
each source carries a serious risk for life of soldiers on a battlefield. 

There is a single condition of BK, which has been viewed as a drawback from the applied standpoint 
for a long time since 1981, This condition is still not lifted. Specifically, BK requires that at least one 
initial condition to be non-zero in the entire domain of interest. However, after getting an extensive recent 
numerical experience with the approximately globally convergent numerical method for MCIPs (see the 
book of Beilina and Klibanov [14] and Section 6 below), the author believes now that this drawback is 
an absolutely insignificant one precisely from the applied standpoint (although the mathematical question 
remains open). Indeed, the most interesting case in applications is the case when the initial condition is 
the function 5 {x — xq) with a fixed position of the source xq. However, replacement of this function by its 
approximation via a narrow Gaussian 5^{x — xq) , 



immediately lets BK working (e > is sufficiently small here). The corresponding boundary data, which 
model the data resulting from a measurement, have only an insignificant change. Therefore, if a numerical 
method for a corresponding MCIP is stable, as it must be, then this change should affect the solution only 
insignificantly. Furthermore, physicists and engineers are indifferent to such a replacement because of the 
above reasons. This is why functions 5 {x — .tq) and 5^ {x — xq) are equivalent precisely from the applied 
standpoint. 

Still, the author has proved [94] uniqueness theorem for an MCIP for the equation uu = Am + 
a (x, y, z) u, (cc, y, z) £ 'E? for the case of an incident plane wave with u (x, y, z, 0) = 0, ut (x, y, z, 0) = 5 (z) 
and under the assumption that this equation for the function u is written in finite differences with respect 
to {x,y) . 

The idea of BK is described in Section 3. Five examples of this section show how BK works. In a less 
general form examples of Sections 3.2, 3.3.2 and 3.4 were first published in originating works [35l[74], also 
see [Ml |82j [86] for more general forms of these three examples as well as for two other examples of Section 
3. In principle, it is possible to formulate BK in a general abstract form, see, e.g. the earlier paper of the 
author [7S] for this form. However, it is not necessary to do so for the understanding of BK. Previously 
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published relevant results are discussed in Section 4 as well as in the end of each of Sections 2,3,5,6. 

2. Carleman Estimates, Holder Stability and the Quasi-Reversibility Method. 

2.1. Definition of the Carleman estimate. We now introduce the notion of the pointwisc Carleman 
estimate for a general Partial Differential Operator of the second order. Let G C M" be a bounded domain 
with a piecewise smooth boundary dG. Let the function ^ € (G) and |V^| 7^ in G. For a number c > 
denote 

{x eG : a^) = c) ,Gc = {x G : i{x) > c) . 
Assume that Gc^ 0- Let Fc C dG a part of the boundary dG defined as 

F, = {x e <9G : <e [x) > c} . 
Let Gc ^ 0. Then the boundary of the domain Gc consists of two parts, 

dGc = diGc U SaGc, diGc = ^„ ^Gc = F^. (2.1) 
Let A > 1 be a parameter, which we will consider to be large. Consider the function tp^ (x) , 

(x) = exp (AC (x)) . (2.2) 

It follows from that 

min ip^ (x) = ip^ {x) |{ = exp (Ac) . (2.3) 

Let A (x, D) be a linear Partial Differential Operator of the second order with real valued coefficients in 
G and with its principal part Aq (x, D) , 

A{x,D)u= aa{x)D'^u, Ao{x,D)u= ^ aa{x)D°'u. (2.4) 

\a\<2 \a\=2 

e G^ (G) for |a| =2,K := max (||a„||pi/^ ^) ; a„ e G (G) for |a| ^ 0,1. (2.5) 

|a|— 2 V \ ^ } / 

Definition 2.1. Let Gc 7^ 0. VFe say t/iat the operator Aq {x,D) admits pointwise Carleman estimate 
in the domain Gc with the Carleman Weight Function ( CWF) ip^^ (x) if there exist constants Aq (i^, K) > 
l,Go = Go {fl, K) > 0,C {fl, K) > depending only on the domain Gc and the number K, such that the 
following a priori estimate holds 

[A^uf ^2 [x)>^Y. {D'^uf pI (x) + GA {Vuf ipl (x) + GA^u^a (^) + div U, (2.6) 

"1=2 

VA > Ao, Vw e G^ (G) , Vx e Gc- (2.7) 
In S2. 6\) the term under the divergence sign satisfies the following estimate 

\U\ < CX' [{\/uf + u'] {x) + ^Y. (^"")' (^) ■ (2-8) 

l"l=2 

In the case of parabolic and elliptic operators Go > and Go = in the case of a hyperbolic operator. 
Lemma 2.1 is elementary. 

Lemma 2.1. Let conditions \2.5]) imposed on coefficients of the operator A be valid. Suppose that the 
Carleman estimate V/i.b]) - ^^^) is valid for the principal part Aq [x, D) of the operator A (x, D) . Then it is 
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also valid for the operator A{x,D), although with a different constant Aq. In other words, the Carleman 
estimate depends only on the principal part of the operator. 
Proof. Wc have 



(2.9) 



where M > is a constant depending only on the maximum of norms Hq-q , |ck| = 0, 1. Substituting 

p.9p in (j2.6|) and taking A sufficiently large, we again obtain (|2.6p . □ 



2.2. Holder stability. Consider the following Cauchy problem for the differential inequality 

\Aou\ < B {\Vu\ + \u\ + I/I) ,Vx e Ge, (2.10) 
I T^^^ 9o{x) ,dnu\T^= gi{x) , (2.11) 



where B ~ const. > and f G L2 (Gc) is a function. Clearly, functions go,9i in (|2.1ip are the Cauchy 
data for the function u. In particular, equation Au = f with the boundary data (|2.1ip can be reduced to 
the problem (|2.10p . (|2.1ip . Wc want to estimate the function u via functions f,go,gi. Such estimates were 
derived in Chapter 4 of the book of Lavrent'ev, Romanov and Shishatskii [112! for parabolic, elliptic and 
hyperbolic operators. 

Theorem 2.1 (Holder stability estimate). Assume that conditions 112. 5\) hold and that the Carleman 
estimate of Definition 2. 1 is valid. Suppose that there exists a sufficiently small number e > such that the 
domain Gc+ss 7^ 0- Denote m = max^ ^ (x) . Define the number (3 = 2s/ (3m + 2e) € (0, 1) . Let functions 
90,91, f be such that go G H^{Tc),9i G L2 (Fc) , / € X2 (Gc) . Let the function u S C^ [Gc) satisfies 
conditions \2.10\) . \2.11\) . Then there exists a sufficiently small number Sq = Sq {e,m, B, K,Gc) G (0, 1) and 
a constant Ci = Ci (e, m, B, K, Gc) > such that if S € (0, Sq), 



II/IIl2(g.) + ll.9o|lHi(r,) + !l.9illL.(r,) < ^ 
then the following Holder stability estimate holds 

<Ci(l + h||^,(e.))'5^V5€ (0,M, 



(2.12) 



(2.13) 



Remark 2.1. Estimate (|2.13p is Holder stability estimate because (3 e (0, 1) . If wc would have /? = 1, 
then (|2.13p would become the Lipschitz stability estimate. Because of the presence of the term ||m||^1(.q ^ , 
this is the so-called "conditional stability estimate" , which are common in the theory of ill-posed problems 
[5l ll41l49ll70[fT45] . Indeed, the presence of this term actually means that we assume an a priori given upper 
bound for the norm |lti||^i(Q j . 

Proof of Theorem 2.1. In this proof G ~ G {e, K, Gc) and Gi = Gi (e, m, B, K, Gc) denote different 
positive constants depending on listed parameters but independent on the function u and the parameter 
A. Obviously Gc+ss C Gc+2£ C Gc+e C Gc. Since Gc+3e 7^ 0, then Gc+2e,Gc+e,Gc ^ 0. Let x[x) be a 
function such that 



xeG^ {Gc),x{x)^ 



Consider the function v. 



1, a; e Gc+2e, 
0,a; G Gc\Gc+e, 

G [0,1] ,2; G Gc+e\Gc+2e- 



XU. 



Then dUni), and I^H^ imply 



\Aov\ < Cl |Vt.| + \v\ + |Vx| |V^.| + I l"l + I/I 

\|a|=2 

^' I = X9o, 9nV |r,= 9odnX + X9i, 
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yx G Gc 



(2.14) 

(2.15) 

(2.16) 
(2.17) 



v{x) ^0,x e Gc\Gc+e- 
Square both sides of (|2.16p . multiply by {x) and apply (|2.6p ignoring Cq. We obtain 



(2.18) 



\a\=2 



^ ' 'uf -divU> 



C\ (^1 - y ^ {S/vf ifl (x) + (^l - wVa (•^) , VA > Ao, V.T e Gc- 
Choose A > Ai := max (Aq, 2Ci) so large that Ci/A < 1/2. Then with a different constant C 



J«l=2 



CA (Vu)^^^ (x) + CA^uVa (2;) ,VA > Ai,V.T € G, 



^ ' -dWU> 



3„,2,„2 



Integrate this inequality over Gc using Gauss-Ostrogradsky formula as well as (|2.ip . (|2.3p . (|2.8p . (|2.14p . 
(PTT)) and dull)- We obtain 



Ga r 
+Ci exp [2A (c + 2e)] 



|Vw|^ +u'^]dx 



(2.19) 



Gc + e\Gc + 2e 



>Ay {Vvf ifildx + J v^ifildx. 

Ga Ga 

Since Gc+3e C Gc+2e C Gc, then strengthening inequality (|2.19l) and using (j2.14p . (|2.15p . we obtain 



Gie^^" J fdx + GiA^e^^™ j \iygof + gl\ dS^ + Ci exp [2A (c + 2e)] 
>A J {\/uf ipldx + J w Vi^a; > A exp [2A (c + 3e)] ^ {\/u f + u^ 
Hence, we have established that 



(\Vu\^ +u'^'j dx 



dx. 



Gie^^™ / fdx + GiX-'e 



3„2A?ri 



(Vffo)' + ff? + Gi exp [2A (c + 2e)] 



> A exp [2A(c + 3£)] 



\H^{Ga + 3e) 



Divide both sides of this inequality by A exp [2 A (c + Be)]. Hence, there exists a number 
A2 = A2 (e, TO, B, K, Gc) > Ai such that 



rdx^ 

[Ga Fa 

Using (|2.12l) . we obtain 



(V.go)' + gl 



dSx 



Gie^^'" + Gl exp [-2Ae] > 



ffl(Gc + 3.) 



2„3A™^^-2Ae||^||2^^^^^ 



,VA>A2. (2.20) 



(2.21) 



We now balance two terms in the right hand side of (|2.2ip via choosing X = X{S) such that 



Hence, 

-2(3m+2e 



A = In (^5""^""+""^ j. (2.22) 

Hence we should have 5 e (0, Sq) , where the number Sq — Sq {e, m, B, K, Gc) is so small that In > 

A2. The target estimate ((2J3)) follows from ^^lT]i and □ 

Theorem 2.2 (uniqueness). Let conditions of Theorem 2.1 hold, in \2.11]) [x) = gi {x) = Q,x £ T ^ 
and also f (x, t) = 0. Then u {x) = for x £ Gc- 

This theorem immediately follows from Theorem 2.1. To prove convergence of the Quasi- Reversibility 
Method (Section 2.5), we need to replace the pointwise inequality (|2.10p with the following integral inequality 

{Aufdx<S^. (2.23) 



Theorem 2.3. Let the function u G [Gc] satisfies inequality li2.2S\) . u \tc= 9„m |rc= and the 
number S e (0,(5) . Assume that conditions 112. 5\) hold and that the Carleman estimate of Definition 2.1 is 
valid. Suppose that there exists a sufficiently small number e > such that the domain Gc+^e 7^ ^- Denote 
m = maxg £.ix). Define the number /3 = 2e/{3m + 2e) G (0,1). Then there exists a sufficiently small 
number Sq = Sq (e, m. A, Gc) G (0, 1) and a constant Gi = Ci (e, to, A, Gc) > such that if S € (0, Sq), then 
the following Holder stability estimate holds 

hllHMG.+3.) < (1 + Mhhg._)) '^^VJ e (0,<5o) • 
Proof. Assume first that the function u E G^ (^c) ■ We have 

^2g2Am > J ^^^^2 ^2 dx > J (Aqu)^ ifl (x) dx - Gi J ((Vu)^ + U^^ ipl {x) dx. 
Gc Gc Gc 

This is equivalent with 

S^e^^"'+CiJ ((Vu)Vm2) ifl {x)dx> j {A^uf ifl (x) dx. 

Gc Gc 



The rest of the proof is similar with the proof of Theorem 2.1. The replacement of m e ipcj with 
u £ [Gc) can be done via density arguments. □ 

2.3. Derivation of the Carleman estimate for a parabolic operator. The goal of this Section is 
to present an example of the derivation of the Carleman estimate. To choose the case of a simplified parabolic 
operator. Our derivation method is similar with the one of §1 of Chapter 4 of the book of Lavrent'ev, Romanov 
and Shishatskii jll2| . 

For any x G M" denote y = (2:2,. .,x„). Let numbers a, r/ e (0,1) and a < ?y. Let F, T > be two 
arbitrary numbers. Consider the function ■!/; (x, t) , 

+ ^ + (2-24) 

Define the domain as 

G^ = {(x, t) : V (a;, i) < ?7, xi > 0} = I xi + M_ + _L_ + a < 7?, 1-1 > 1 . 



Let A, > 1 be two large parameters which we will choose later. Consider the function Lp^^ ^ (2;,^) , 

V9(a;,t) = exp(AV^-'') . (2.25) 

Lp^^ ^ {x,t) is the CWF for our parabolic operator L (below). To simplify notations, we use the notation 
ip{x,t) instead of ip)^ ix,t) . Hence, the boundary of the domain G,, consists of a piece of the hyperplane 
{xi = 0} and a piece of the paraboloid {ip {x, t) = cci > 0} , 



|2 ,2 



Consider a function a {x,t) for {x,t) G G,, such that 

a S (Gri) , K = 11111(71(0 ) ; (2^1 ^) 5^ '^o = const. > for {x, t) S G^. 
Consider the parabolic operator L, 

Lu ~ ut — a (x, t) Au, (x, t) € G,j. 



(2.26) 
(2.27) 

(2.28) 
(2.29) 



By Definition 2.1 we want to estimate now {Luf ip"^ from the below. Introduce the new function v = wp 
and express derivatives of the function u via derivatives of the function w, using p.24p and ()2.25p . Below 
O (1/A) , O (1/j/) denote different G^ {G7]) —functions, which are independent on the function u, and such 
that |0(1/A)| < C/A, |0(1/^)| < Gjv.'iX.v > 1. Here and below in this section G = C{ao,K,G,^) denotes 
different positive constants depending only on listed parameters. 

Lemma 2.1. Suppose that the function a{x,t) satisfies conditions 12. 28]) . Then there exist sufficiently 
large numbers Aq = Aq (cq, Gn) > 1, J^o = (oq, K, Gn) > 2 such that for any function u G G^'^ (Gjj) the 
following estimate holds for all A > Aq,;^ > j/q, {x,t) G G^ 



3,,4„;,-2iy-2^2 



{Luf ij^+^ip^ > ^GXv {Vuf ip^ + GA^z/V 
\Ui\ + \Vi\< GA^ 1.3^.-2.-2 ^^^^^2 ^^2^ ^; 

Proof. We have u = v exp (— AV'"'^) ■ Hence, 

(vt + -^Xvip'^'^v^ exp (-AV'"") , 
[vx^ + Ai/V'"''-^'y) exp (-AV'-'^) , 



div Ui + 9tVi, 



(2.30) 
(2.31) 



1 + 



exp (— At/i 



(^^i + ^^^^ ^u^exp(-AV' ''),ie[2,n] 



These equalities imply that 



a"" 



exp (-A?A , i e [2, n] . 



Ir ^2 ,1^+2 2 

(Lu) ip ip 



I't — aAv 



2aXv'ip ^ "^Vx 



i=2 



Y2' 



2,,2„/,-2iy-2 



aA^i.> 



E 

1=2 



y4 



1/+2 



Denote 



zi = Vt,Z2 = ~aAv, 



Z3 



Z4 = —aX v 4' 



2,,2„/,-2!y-2 



1 + E 



n 9 



i=2 
1 



i=2 



Hence, 



{Luf ilj^+^ifi^ > {zf + 2ziZ2 + 2ziZ3 + Z3 + 2z2Z3 + 2Z3Z4 + 2Z1Z4) 7/-''+^ 
Step 1. Estimate 22122-0"^^ from the below. By ([QS]) 



u+2\ 



1=1 



+2 ^ V,, a0"+' + 2wt ^ v^^ (aV'"+') 



i=l 



1=1 



1.^+2 



{Vvf {ai'''+^)^ + 2vtJ2^-. H' 



i/+2\ 



i=l 1=1 

Thus, since < I7 then 

n 

2^1^20"+' > 2z;t - (^«)' + t^i,i + (^i.i)t : 

j=i 

Step 2. Estimate (zj + 2ziZ2 + 2ziZ3 + z|) i)"^'^ . Using ((O^ . ([OS)) and (P37| . we obtain 

(z^ + 2ziZ2 + 2ziZ3 + zl) > 



4 + zl + ' 



- 2zi j^Z3 + ^ v^^ j - C^i' (V«)' + div Ui.i + (Fi.i), 

>zl + zl-zl- |^^3 + ^^^., (a0"+'),^^ - Ci.(Vi;)' + div C/i,i + (1/1,1), 

n / n \ 2 

= - - 2z3^^;,. ("V'"^'"').. - ( E^-^ («^''^')., I - Cz/(Vi;)' +divC/i.i + (^1,1 
1=1 

Using (|2.33p . we obtain 

n 

-2z3 Y («^''^').. = ^a^Ai/ (jy + 2) + Y Tk^x 



,4=1 



i=l 



i=2 



- CAi/ (Vw)^ > -CXv (Vw)^ . 



Thus, 



{zl + 2ziZ2 + 2ziZ3 + zl) > -CAj/ (Vi;)^ + div Ui.i + (Vi,i), 



step 3. Estimate 2z2^3'0'' > 
Consider for example the term Aa^Xvipyx-^^/^v, 

n n 

i=l i=l 
n n 

- ^ AXiya'^ipiyxiVxiXi - ^ 4Ai/ (a^V')^. i^xi'yxi 



(2.39) 



1=1 



i=l 



> {AXva^Vx,Vx,i^)^, + ( ^ 2Ai/aVj^I, ) - '2Xv (aV)^^ (Vw)^ - CAi^ {'sjvf 

1=1 \i=l / a;^ 

> -CXv {\Ivf + AwUi.2. 

Thus, 'ia'^XvipVx-i^^v > —CXv (Vu)^ + divC/i.2. Similarly we obtain using (j2.39p 

222^31^''+^ > -CXv {Vvf + div C/i,3. 



(2.40) 



Step 4. Estimate 2z3Z4-0'^'^^, 

o ; 1^+2 A 2 \ 3 3 / -2iy-l 



"2 /I 



1=2 



\ i— 2 / J xi 



n 2 



i=2 



y4 



E 

^j=2 



y2 



+2a^AV'^ (2j/+ 1) V" 



1+ El 



y2 



o 



o 



Thus, 



Similarly, 



22324-0"+^ > CA^i^V"^''"''^^ + divi7i,4 



3,,4„;,~2i/-2„,2 



2ziZ4^''+^ > -CA^j/^^/'^'^^'w" +divC/i,5 + (V^i,2)t 



, 2,,3„/,-i^-l„,2 



(2.41) 



(2.42) 



Summing up (|2.38p - p.42p and replacing v with u = vLp~^ , we obtain p.30p . (|2.3ip . □ 

As one can see that we have one positive and one negative term in the right hand side in estimate p.30p . 
Therefore, we need to balance them somehow to obtain only positive terms. To do this, we prove Lemma 
2.2 first. 

Lemma 2.2. Suppose that the function a{x,t) satisfies conditions \2.28\] . Then there exist sufficiently 
large numbers Aq = Aq (oq, K, Grj) > IjVq — vq (oq, K, G,,) > 2 such that for any function u G C^'^ (^'?) 
following estimate holds for all X > Xq , v > h'o i [x , t) (z Gn 



{ut - a (x, t) Au) u(f'^ > C {Wuf Lp^ - CA^z/^V"''''~V^w'' + div U2 + dtV^ 
\U2\ + \V2\ <c(|Vu|'+u') 

10 



(2.43) 
(2.44) 



Proof. We have 

{ut - a (x, t) Au) = i^^^'^j -^^'^'^'"'^f'^'^'^ + X] ^) + i^^f '^'^ 

^ ' * j=l 

n n 

—2a\v 'ip^.'ip~^~^UxiU(p^ + ttxiUxiUip^ 

1=1 4=1 

> C (Vu)^ vs^ - CA^^^V^'^^'V^^t^ + div f/s + V^2t- □ 

Theorem 2.4. Suppose that the function a {x,t) satisfies conditions \2. 28]) . Then there exist sufficiently 
large numbers Aq = Aq (ag, Gn) > 1, i^o = ^^o (ctOi -^'i G,,) > 2 smc/i t/iat for any function u G C^'^ (^'j) 
following Carleman estimate holds for all A > Ao,i^ > i^oi (^^lO G 

(Lu)^ > (jxi, (Vu)^ _|_ C'A^z/V~^''~V^w^ + div [/ + Ft, (2.45) 
|C/| + < CA^J^V^"""" ((Vu)' + u^) ip^. (2.46) 

Proof. Multiply ((2:43| and ((2:44)) by 2Ai/ and sum up with ([230)) . ((23T|) . We obtain 

(iw)^ i/'^^V^ + 2CAi/ (Lu) > CAi/ {\7uf Lp^ + CA^i^V'"^""^ (^1 - ^ f^^i^ + div [/ + Ft, (2.47) 

where the vector function {U,V) satisfies (|2.46p . Choose = {ao,K,G,f) > 2. Also, since ip'^^'^ < f, we 
have 

{Luf ij^+^ip^ + 2\v (Lu) u(p^ < 2 {Luf (p^ + X^v'^p^u^. 

Combining this with (|2.47p . we obtain p.45p . □ 

We now want to incorporate higher order derivatives ut, u^^x in the Carleman estimate of Theorem 2.3. 
To do this, we prove Lemma 2.4 and Theorem 2.4. 

Lemma 2.4. Suppose that the function a{x,t) satisfies conditions i2.28\) . Fix the number v :~ 
1^0 {ao, K,Gri) > 2 of Theorem 2.3. There exists sufficiently large number Xq = Xq {oq, K,Gri) > 1 such 
that for any function u £ i^v) following estimate holds for all A > Aq, {x, t) € G,, 



- CA^ ( Vw) + div [/3 + 9*1^3, (2.48) 



|C^3| + \V3\ <c\^u',+ ul,^ + (Vufj p\ (2.49) 
Proof. We have 

{Luf p^ > [u1 + 2autAu + {Auf^ p^ . (2.50) 
Step 1. Estimate 2autAu(p^, 

n n n 

2autAup'^ = 2autUxiXiV'^ = (2awtu^;(y9^)_^ — 2 autxjUxi'p'^ 

i—l i—1 i—1 

n n 

-2'YaxiUtUx,(p'^ - 4Ai''0"''""' i4')x^ autUx,(p'^. 
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We have 

Since the number v :~ vq (ao, K, G,,) depends on the same parameters as the constant C, we can incorporate 
V in C . Hence, we obtain 

2aut/^ULp^ > -u^if^ - CX^ {Vuf ip^ + div U^.,! + (Vg.O^ . (2.51) 
Step 2. Estimate (Aw)^(p2, 



n / n 



*J = 1 j = l \»=1 / Xj 

n n n 
-^^Ux,XiX,Uxjip'^ +2Xvip~''~^ ^ UxiX,Ux^^XjV>^ 
i=l j = l i,j=l 
n I n \ \ ^ 

^=1 V ^ = 1 / x^ = l 

n / n \ 1 

^ ^w..x.wx,<p' = 2 ^ <,^(^2-CA2(VM)'^2^divC/3,2, 

1 \i=l / Xj hJ = l 



2 2 



Hence, we have obtained that 



1 " 

(A«)^ V^'>^T. (^")' + ^3,2, (2.52) 

■i j' — 1 

|C/3,2| ^ ^ E + (2-53) 

Comparing JlSOl), ([23T|) . ([232| and ([233]) . we obtain ^(T^ and (|2^ . □ 

Theorem 2.5. Suppose that the function a{x,t) satisfies conditions 12.28\) . Fix the number v 
j^o (flOi A, G,,) > 2 of Theorem 2.3. Then there exists a sufficiently large number Aq = Aq (oq, isT, G,,) > 1 
such that for any function u ^ ipri) the following estimate holds for all A > Aq, [x^t) S Gn 

{Luf ip^ >^{ul+Y^ ul^^^ I (^2 ^ GA {\Iuf ip^ + GA^uV^ + div C/ + Vi, (2.54) 



A 

Co 



\U\ + \V\<'^ \u^, + J2 + CX' {{^uf + u') p\ (2.55) 

where the constant Gq = Gq {oq, K,Gri) > depends on the same parameters as the constant C above. 

Proof. Divide incquahty p.48p by 6 A with an appropriate constant h ~ b {uq, K,G,f) > 1. Next, add 
the resulting incquahty to (j2.45p . where set v i^o- Then we obtain (j2.54p . p.55p . □ 

Remark 2.2. If one would set above ut = and would ignore the term /T'^ in the function ijj, then 
one would obtain analogs of these results for the general elliptic operator. A close analog of Theorem 2.4 
is valid for a general parabolic operator of the second order, see §1 of Chapter 4 of the book of Lavrent'ev, 
Romanov and Shishatskii [112) . Therefore, Theorems 2.1-2.3 are valid for this operator. 
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2.4. Carleman estimate for a hyperbolic operator. Theorem 2.5 is proven in tlie book; of Beilina 
and Klibanov [13] . Similar theorems were established in earlier books of Isakov [S51 [57] and Klibanov and 
Timonov [S5]. Theorem 2.6 can be found in §4 of Chapter 4 of the book of Lavrent'cv. Romanov and 
Shishatskii |112j . We do not reproduce proofs here for brevity. 

For brevity we consider a simple domain = {|a;| < R} C R". Let T = const. > 0. Denote 

Qt = X (0, T) , Q± = r2 X (-T, T},ST^dnx (0, T),S$ ^d^x (-T, T) . 

Choose a point xq £ R". In particular, we can have Xq G 17. Let the number 77 e (0, 1) . Let A > 1 be a large 
parameter. Define functions ^ (x, t) , ip {x, t) as 

^{x,t) = \x-xo\^ -vt^,ipix,t) = exp [A^ (x, t)] . (2.56) 

For a number 7 > define the domain Gj as 

G, ^ {ax,t) > j} . (2.57) 

One can choose 7 such that 

G-y C Qt- (2.58) 

Theorem 2.5. Lei n = < i?} C R",n > 2. Lei conditions f05|) . ifg^PjJ and fHig) 6e m pZace. 
Lei 0? = const. > 1. Lei i/ie function c{x) satisfies the following conditions 

(x) e ,Vx e n,c e (H) , (2.59) 
(x- a:o,Vc"^ (x)) > 0, Vx e n, (2.60) 

/or a certain point xq G f2, where (•,•) denotes the scalar product in R". Lei Lm = utt — (x) Am be a 
hyperbolic operator. Let 

P = P(xo,fi) = max|x- xol . (2.61) 

Then there exists a number r/Q = r/Q ^17, d, P, |j Vc||p^j^-j^ G (Ojl) such that for any rj G (0,7]g) one can choose 
a sufficiently large number Xq = Xq (^d, P, || Vcjlpj-jj-j , 77, 7^ > 1 and the number G ~ G (^d, P, II Vc||pj-j^-j , 77, 7^ > 
0, such that for all u £ G'^ ifii) '^'^^ /''^ A > Aq the following pointwise Carleman estimate holds 

{Lufip^ > GX (\Vu\^ +ufj (p^ + X^u^ip^ + diYU + Vt, in G^, 

where 

|C/| < CA^ (IVmIVu^ + u^) (2.62) 
|y| < CA^ [|i| (uj + |Vn|^ + u^^ + (|Vn| + |u|) |Mt|] (2.63) 

In particular, \2. 63\) implies that if either u (x, 0) ~ or Ut (x, 0) = 0, then 

y(x,0) = 0. (2.64) 



Theorem 2.6. Let c (x) — 1. Then Theorem 2.5 is valid for any rj € (0,1) 
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2.5. The Quasi-Reversibility Method (QRM). In this section we use notations of Sections 2.1, 
2.2. Let A be the operator of Section 2.1. QRM delivers an approximate solution of the following Cauchy 
problem 

Au = f,xeGc, (2.65) 
/ e L2 {G,),ue Hi {Gc) ^ {u e (G^) : u |r,= d^u |r,= O} . (2.66) 

To find that approximate solution, QRM minimizes the following Tikhonov functional with the regularization 
parameter 7 

(u) = \\Au - /lli^(e^) + 7 \\u\\l.^a.^ ,ueHl (Ge) . (2.67) 

The variational principle implies that any minimizer G Hq (Gc) satisfies the following integral identity 

iAu^,Av) + 7 [u, v] = (/, Av) yveH'l (G^) , (2.68) 

where (, ) and [, ] are scalar products in L2 (Gc) and (Gc) respectively. Riesz theorem and p.68p imply 
Lemma 2.5. 

Lemma 2.5. Let A be the operator defined in i2.4\ l, 12.5]) and the function f € L2{Gc)- Then for 
any 7 > there exists unique minimizer u £ (Gc) of the functional 1^2.67^ . Furthermore, with a constant 
C = C (A) the following estimate holds \\u-y\\jf2(^Q j < 

Theorem 2.7 (convergence). Assume that conditions 12. 5\) hold and that the Carleman estimate of 
Definition 2.1 is valid. Suppose that there exists a sufficiently small number e > such that the domain 
Gc+3e 7^ 0- Let the function u* be the exact solution of the problem i2.65\) . 112. 66\) with the exact function 
f* G ^2 (Gc) . Let 11/ — 11^2(0 ) — ^ '^'^'^ ^ ^ ■ Denote m = max^ (x) . Define the number j3 ~ 
2e/ (3?7i + 2e) G (0,1). There exists a sufficiently small number Sq = Sq ^e, m, A, Gc, ||m*||^2(q G (0,1) 
and a constant Ci = Ci (e, m. A, Gc) > such that if S d (0, Sq), then the following convergence rate is valid 



<Ci{i + \\u*\\^,^^^^^s^/^yse{o,So). 



Proof. We have 



{Au* ,Av)+^[u*,v] = {f*,Av)+-f[u*,v],yveHl{Gc). 
Subtract this identity from (|2.68p and denote ~ Uj ~ u* , g ^ f — f* . We obtain 

iAw^,Av)+-,[w^,v] = {g,Av)~-,[u*,v],yve HliGc). (2.69) 
Setting in (j2.69p v :~ and using Cauchy- Bunyakovsky inequality, we obtain 

ll^«^7llL(Ge) - + ll^*llH^(Ge)) ^ 

||w7ll^2(G^) < i + I|w*|Ih2(g.)- 

The rest of the proof follows from Theorem 2.3. □ 

2.6. PubUshed results about QRM. As it is clear from (|2.67|) . QRM is a special form of the Tikhonov 
regularization functional, see [5j [Ml |49l I145j for the theory of this functional. In "conventional" Tikhonov 
functional the originating operator is continuous. On the other hand in QRM the originating operator is a 
PDE operator, which is continuous only if its domain is H^. QRM is well suitable for providing approximate 
solutions for ill-posed Cauchy problems for PDEs, including boundary value problems with over-determined 
boundary conditions. QRM was first introduced by Lattes and Lions in their book |111| . This book shows 
how to apply QRM to ill-posed Cauchy problems for all three main types of PDE operators of the second 
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order: elliptic, parabolic and hyperbolic. Although convergence theorems were proven in convergence 
rates for QRM were not established there. 

The first work where a Carleman estimate was applied to get convergence rate of QRM was one of 
Klibanov and Santosa [80]. In [80] QRM was applied to the Cauchy problem for the Laplace equation. In 
Chapter 2 of the book of Klibanov and Timonov [86] Carleman estimates were also applied to establish 
convergence rate of QRM for ill-posed Cauchy problems for elliptic, parabolic and hyperbolic PDEs. Next, 
Carleman estimates were used to prove convergence rate of QRM for the Cauchy problem for the Laplace 
equation by Bourgeois [32], Bourgeois and Darde [33] and Cao, Klibanov and Pereverzev [IT]. The QRM for 
the problem of determining of the initial condition in the parabolic PDE from boundary measurements was 
considered by the author in [89] [95]. Papers [33l|4Tl|80] contain numerical results. As to the application of 
QRM to the problem with the lateral Cauchy data for the hyperbolic PDE, see Sections 5.4 and 5.5. 

While above citations of QRM are concerned only with linear problems, it was recently applied by the 
author with coauthors to solve MCIPs with backscattering data via the approximately globally convergent 
method [Ml HOOl 11011 11021 1103] , also see chapter 6 of [H] . The main difference between the latter application 
of the QRM and the conventional one is that MCIPs are nonlinear. 

3. The Bukhgeim-Klibanov Method. 

3.1. Estimating an integral. First, we estimate a Volterra-like integral with a weight function. For 
the first time an analog of Lemma 3.1 was proven by the author in [73]. Next, that proof was published in 
some of above cited follow up papers of the author about BK. The estimate of this lemma with the parameter 
1/A in it was first published in the book of Klibanov and Timonov [86], also see Section 1.10.3 in the book 

M- 

Lemma 3.1. Let the function ip ^ [0,a] and (p' (t) < —b in [0,a], where b = const > 0. For a 
function g €z L2 (—a, a) consider the integral 

I{g,X)^J (^J^g{T)dT^ exp [2Xip{t'^)]dt,X^ const. >0. 



The 



/ (5, A) < ^ y 5' it) cxp [2A^ {t^)] dt. 



Proof. We have for t > 



texp [2X^{t')] ^ii^^exp [2X^ {t')] 



' {exp[2A^(t^)]} = -— ^4{-exp[2A^(i^)]} 



Hence, 



AXp>'{t^)dt ^ yjj AXip'{t^)dt 

f g (r) dr exp {2Xif (t^)) dt< exp (2A^ (t^)) t (^J^ g^ (r) dr^ dt 
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cxp (2A^ {a')) £ (r) dr + ^ £ (r) cxp {2X^ (t^)) dt 
g'^ (T)exp {2Xip {t"^)) dt. 



1 

- 4A6 

Thus, we have proved that 
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exp(2A^(t2)) (^^ 5 (r)dT^ dt < ^ g^r) cxp {2X^ (t^)) dt. 



Similarly, 



cxp {2X(p (t^)) (^J^ g (r) dr^ dt < ^ j g^ (r) exp {2Xip (t^)) dt. □ 

3.2. An MCIP for a hyperbolic equation. In this section we use notations of Section 2.4. Let 
functions Ua {x,t) E C {Qt) i |q^| < 1 and the function c{x) G (fj) ,c(a;) > const > 0. Let the function 
u & [Qrp) be the solution of the following initial boundary value problem 

c(x)utt = Au+ ^ aa{x)D'^u, in Qt, (3.1) 

|al<l 

u{x,0) = fo{x), ut(x,0) = /i(x), (3.2) 
u\s^^p{x,t),ST^dnx{0,T). (3.3) 

Coefficient Inverse Problem for the Hyperbolic Equation (|3.ip . Let the Neumann boundary 
condition be known, 

— \sr^qix,t). (3.4) 

Determine one of x— dependent coefficients of equation (jSHJ, assuming that other coefficients are known, so 
as and functions fo,fi,p,q in p.2p - p.4|) . 

Theorem 3.1. Let the domain Q = {\x\ < R} C > 2. Denote b{x) = Xj \fc(x). Let the function 

b (x) satisfies conditions 112. 5 9\) . 112. 60\) . In addition, let coefficients Oa G C (fl) . Consider two cases: 

Case 1. The coefficient c{x) is unknown and all other coefficients aa {x,t) are known. In this case we 
assume that 

A/o (x) + 5Z (^) ^"/o (^) ^Oforxen. (3.5) 

|q|<1 

Then for a sufficiently large T > there exists at most one pair of functions {u, c) satisfying iS. 1\) - (K4^ and 
such that u E {Qt) ■ 

Case 2. Let ao be one of multi-indices in h3.1\) . |q;o| < 1- Let the coefficient aa„ {x) be unknown and all 
other coefficients are known. In this case we assume that 

D^''fo{x)^OforxeTi. (3.6) 



Then for a sufficiently large T > there exists at most one pair of functions (u,aao) satisfying ll3.1\) - (3^ 
ana such that u e C^+l""! (Q^). 

// in is. 2\) fo (x) = 0, then conditions of these two cases should be imposed on the function fi [x), the 
required smoothness of the function u should be u ^ {Qt) Case 1 and u g C4+I"«l {Qj) in Case 2, 
and the above statements about uniqueness would still hold. 
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Proof. First, we note that if /q (x) = 0, then one should consider in this proof Uf instead of u, and 
the rest of the proof is the same as the one below. We prove this theorem only for Case 1, since Case 2 is 
similar. Assume that there exist two solutions (ui,ci) and (1*2,02). Denote u = ui — 1*2, c = ci — C2. Since 



ClUltt 

then KB-^M lead to 



C2U2tt = ClUltt - ClU2tt + (Cl - C2) U2tt = CiUtt + CU2tt, 



Lu ~ Cl {x) utt — Alt — {x) D"u = c{x) B {x, t) , in Qt, 

u{x,0) = 0,ut {x,0) = 0, 
~i '^"1 n 

u\St = ^Ist = 0, 

B {x,t) : = -U2tt {x,t) . 
Setting in c '.^ C2,u U2,t :— and using (|3.5p and (|3.9p . we obtain 



(3.7) 

(3.8) 
(3.9) 
(3.10) 



B (x, 0) = -C2' (x) A/o (x) + ^ flc. (x) DZh for x e f!. 

V i"i<i 

Hence, there exists a sufficiently small positive number e, such that 

B(a;,t) 7^ in =IIx [0,£]. 

By dUl) 



Hence, 



Or 



Denote 



By ([XTT]) and (jXTOl) 



Denote 



c{x) 



Lu 



B{x,t) 



in 



d d 



B{x,t) 



= in 



Lut = (Lu) in 



^(a;,i) = (a;,i) ■ 



/i e C2 (Q^ 



V {x, t) — Ut (x, t) — /iw (x, t) . 



(3.11) 



(3.12) 
(3.13) 

(3.14) 
(3.15) 



We can consider p.lSp as an ordinary differential equation with respect to u [x, t) with the initial condition 
from (??), i.e. u{x,Q) = 0. Hence, using (|3.13p . (|3.14l) and (|3.15p . we obtain 



u{x,t)=^ J K {x,t,T) V {x,t) dr, 


K{x,t,T) = eC^Tix [0,e] X [0,£]) 



Bix,T) 



v{x,0) = 0. 



(3.16) 

(3.17) 
(3.18) 
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Using (IS331), (|XTi)) . (|XT5|) . (pUj) and (PT7|) . we obtain in 



Ci (ut)j( - /iCiUit = civtt + 2cihtUt + cihttu 



t t 
ciVtt + 2cihtv + 2ciht J Kt {x,t,T) v {x,t) dr + cihtt J K {x,t,T) v {x,t) dr, 





Aut - hAu = Av + 2VhVu + Ahu 

t 



Since by (|3.12p and p.l3p Lut — h ■ Lu — in Q^, then two recent formulas, boundary conditions p.Sp . the 
initial condition p.lSp as well as p.l4p and p.l7p lead to 



|ci (x) vtt - Av\ < M 



\Vv\ {x,t) + \v\ {x,t) + J i\Vv\ 



dv 



\v\) (x,t) dr 



V I 5. = ^k=0, 



v{x,0) = 0, 



(3.19) 

(3.20) 
(3.21) 



where M > is a constant independent on v,x,t. The idea now is to apply the Carleman estimate of 
Theorem 2.5 to the problem p.l9p - p.2ip and estimate integrals using Lemma 3.1. 



Let the point xq € ri, the number P ~ P (xq, fi) be the one defined in p.6ip and ijq = ijq yd, P, || Vc|lp^j^-j j G 
(0, 1) be the number considered in Theorem 2.5. Choose an arbitrary number rj G (O,77o) . Assuming that e 
is so small that 770^^ < consider the domain G^^s, 



Then G^^2 C Q^. Let Se ~ dfl x (0,£) . The boundary of G^^2 consists of three parts, 



a2G+, = {|x|e (Vi?2-,ye^i?),t = o}, 
^3G+2 = {|x|' - 7jt^ = i?2 - Tje^,t > 0, |x| < i?} 



Square both sides of inequality p.l9p . multiply by the function ip'^ {x,t) defined in subsection 2.4, apply 
Theorem 2.5 and Gauss-Ostrogradsky formula. By (|2.62p and p.20p integral over diG~^ 2 equals zero. By 
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(|2.64p and (|3.2ip integral over d2G^ 2 also equals zero. Hence, we obtain with a different constant M 



C\ / {^v) Lp^dxdt + Cy / v-'tp-'dxdt < M / (Vv) +v 



(p^dxdt 



G+2 



-CX^ exp [2A - Tje^)] I {Vvf + 



M 



dSx,t,yx > Xq. 



J y"(|Vf| + b|)(a;,r)dr ip'^dx^.22) 



For an arbitrary point (x, 0) G d^G^^ 2 consider the straight line, which is parallel to the t— axis and 
passes through the point (a;,0). Then this line intersects the hypersurface ^sG"*" 2 at the point {x,t{x)). 

Hence, applying Lemma 3.1 to any function g (z C (^G^^2^ , we obtain 

J i J \g\{x,T)dT\ (f>'^dxdt= J dx J i J \g\{x,T) dr] ip^dt <J J g'^V^dxdt. (3.23) 



^+ \o 



03 



\0 



Using p.22p and (|3.23p and choosing a sufficiently large number Ai = Ai (A/, Ao) > Aq, we obtain with 
a different constant C 



CX J {VvYif^dxdt + GX^ j v'^ip^dxdt<GX^Gii^[2X{R^ --qe^)] 

G+ 2 G+ 2 d3G+ 



Let 5 e (0, ?ye^) be any number. Then C G^^^2 and ip'^ {x, t) > exp [2A (i?^ — 6)] in G^. Hence, strength- 
ening the last inequality, we obtain with a different constant C 



X^exp[2X{R^ -S)] J v^dxdt<GX^cxp[2X{R^ -rie^)] J {\/v f+v 
Dividing by A^ exp [2 A (i?^ _ S)] , we obtain 



dSxf,yX > Ai. 



93 G+ 



v'^dxdt < C exp [-2A (r/e^ - 5)] / (Vi;)^ + 



dSx, 



Setting here A — 00, we obtain v (x, t) = in G^. Since numbers 77 G (0, ?7q) and S G (0, 776^) are arbitrary 
ones, then v (x, t) = in G 



Hence, by p.l6|) u {x, = in G^ ^2 • Hence, setting i = in p.7p and using p.Sp . we obtain 



■(x) = for x e £ (/R^— 



(3.24) 
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Substitute (jX^ in dXT]) and use (??) and dS^). We obtain for {x, G | e (^y/ - rj^e^ , i?) | x (0, T) 

n 

Lu = ci (x) utt — Am — ^ Ua (x) D^u — 0, (3.25) 
i=i 

u(a;,0) =0,ut(x,0) =0, (3.26) 
u\sr = -g^\sT=0- (3-27) 
Consider an arbitrary number to G (0,r — e) . And consider the domain Gij^ea (ip) , 

Gvo^^ (to) = {{x,t) : \x\^~iio {t-tof > R^ -T]o£^,t>0,\x\ < i?} . 
Since to £ {0,T — e) , then t E (0, T) in this domain. Hence, 

G,^,2 (^o) C { e (/R^^^^, i?) } x (0, T) . 

Hence, using conditions p.26p . p.27p . we can apply Theorems 2.3, 2.5 to the domain G^^^2 (to)- Therefore, 
u{x,t) = for {x,t) € C!,j^^2 (to) ■ Since to is an arbitrary number of the interval (0,r — e) , then, varying 
this number, we obtain that 

u{x,t) =0 for {x,t) e G (^^/R^~^,R^^ x {0,T - s) . 

Therefore, we now can replace in p.7p -(??) sets Qt and St with sets Qj^ and Sj^ respectively, where 

Q't = {\^\ < /R^"^^} X (0,r - £) , = = ^R^-r^oE^} x (0,T - e) . 

Next, we repeat the above proof. Hence, we obtain instead of (|3.24p 

c(x) = for X e e (/R2r^2^,i?)| . 

Since e > is sufficiently small, we can choose e such that R^ = krjoe^ where k = k (R, e) > f is an integer. 
Suppose that 

T>ke^—. 

Then we can repeat this process k times until the entire domain = < R} will be covered. Thus, we 
obtain that c{x) = in 17. Hence, the right hand side of equation p.7p is identical zero. This and the 
standard energy estimate imply that u{x,t) = in Qt- D 

An inconvenient point of Theorem 3.1 is that the observation time T is assumed to be sufRciently large. 
An experience of the author of working with experimental data |141 ITSl f92l I102| I103j indicates that this is 
not a severe restriction in applications. Indeed, usually the pre-processing procedure of the measured signal 
leaves only a small portion of the time dependent curve to work with. Still, it is possible to restrict the value 
of T via imposing the condition fx [x) = 0. The proof of Theorem 3.2 partially uses arguments of works of 
Imanuvilov and Yamamoto [601 161) . 

Theorem 3.2. Assume that all conditions of Theorem 3.1 are satisfied. In addition, assume that the 
function fi (a:) = 0. Then Theorem 3.1 holds if 



(3.28) 



In particular if c{x) = 1, then it is sufficient to have T > R. 

Proof. Just as in the proof of Theorem 3.1, we consider now only Case 1, keep notations of Theorem 
3.1. Introduce the function w {x,t) = uu {x,t) . Then by (|3.7p -(??) 



Cl (x) Wtt - Aw - ^ fla (x) D'^w = -c (x) dfu2., in Qt-, (3.29) 



(3.30) 



By (p:^ and ([53T|) 



-c{x) 



wt (x,0) = 


0, 




dw 1 


w\St = 


dn ^ 


w {x, 0) = 


-c(x)j 


p{x) = 


cr^ (x) 


w (x, 0) 


1 


p(x) 


p{x) 



(x) I A/o (x) + ^ (x) i?^/o (x) I =^ 0, X e r! (3.31) 

a|<l 



t 

w{x,t)— / Wt{x,T)dT 







Substituting this formula in p.29p and proceeding similarly with the proof of Theorem 3.1, we obtain that 
c (x) = in and w (x, t) ~u (x, t) = in Qy. □ 

3.3. MCIPs for parabolic equations. In this subsection we prove uniqueness theorems for three 
MCIPs for parabolic PDEs. Unlike the hyperbolic case, conditions imposed on the principal part of the 
corresponding elliptic operator are not restrictive in two out of three of these problems. The reason is that 
Carleman estimate can be proven for a quite general parabolic operator, see §1 of Chapter 4 of the book 
[112) . Below C'fe+zs^ C'2fe+,3,fe+,3/2 ^^^^ Holder spaces, where fc > is an integer and /3 G (0, 1). 

3.3.1. The first MCIP for a parabolic equation. Denote Dlf+^ = M" x (0, T) . Consider the Cauchy 
problem for the following parabolic equation 

c (x) ut^ An + aa,{x)D'^u in Dp-\ (3.32) 

|q|<1 

«(x,0) = /o(x), (3.33) 
c, a„ e C" (M") , c (x) e [1, d] , foe C^+^ (R") . (3.34) 

The problem p.32p - p.34p has unique solution u e (^2+^,i+^i/2 (^Jf^^^ ^ g^e the book of Ladyzhenskaya, 

Solonnikov and Uralceva [104j . Just as in Section 3.2, we assume that O = {|x| < i?} C IR",n > 2. Let 
r C be a part of the boundary of the domain Vl,T = const. > and Tt = F x (0, T). 

The First Parabolic Coefficient Inverse Problem. Suppose that one of coefficients in equation 
(|3.32p is unknown inside of the domain f2 and is known outside of it. Also, assume that all other coefficients 
in (|3.32p are known and conditions p.33p . p.34p are satisfied. Determine that unknown coefficient inside of 
f2, assuming that the following functions p (x, t) and q (x, t) are known 

Ou 

u|r^=p(x,0, —\rT,= qix,t). (3.35) 

It is yet unclear how to prove a uniqueness theorem for this CIP "straightforwardly" . The reason is that 
one cannot extend properly the solution of the problem p.32p . p.33p in {t < 0} . Thus, the idea here is to 
consider an associated MCIP for the hyperbolic PDE p.ip using a connection between these two CIPs via 
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an analog of the Laplace transform. Next, since the Laplace transform is one-to-one, then Theorem 3.1 will 
provide the desired uniqueness result. 

That associated hyperbolic Cauchy problem is 



c (x) vtt = Aw + ^ flc, (x) Dy in DJJ^+i = M" x (0, oo) , (3.36) 

|q|<1 

u|t=o = 0, vt\t=Q = /o [x) . (3.37) 

In addition to (|3.34p , we assume that the coefficients c {x) , {x) and the initial condition /o (x) are so 
smooth that the solution v of the problem (|3.36p . p.37p is such that (a) v € \ D^] if the function 



c{x) is unknown, and (b) v G C^+I"l yD^j if the function c{x) is known and any of functions {x) is 
unknown. 

Consider the Laplace- like transform 72, of Reznickaya [132j . Since the publication jl32) in 1973 this 
transform is widely used |14[ 1861 1112] . The following connection between solutions u and v of parabolic and 
hyperbolic Cauchy problems p.32p . p.33p and p.36p . (|3.37p can be easily verified 



u [x, t) 



2Wnt 



poo 




/ exp 

Jo 


4t_ 



Tv{x,T)dT ■.= nv. (3.38) 



Using an analogy with the Laplace transform, one can easily prove that the operator TZ is one-to-one. Hence, 
given functions p, q in p.35p . the following two functions p{x,t) and q{x,t) can be uniquely determined 

_ dv _ 

v\r^=p{x,t), — |r^=g(x,i). (3.39) 

Therefore, we have reduced the First Parabolic Coefficient Inverse Problem to the hyperbolic CIP p.36p . 
(13371) . (13:39)) . 

To apply Theorem 3.1, we need to replace Foo hi p.39p with Soo- To do this, we observe that, using 
p.35p and the fact that the unknown coefficient is given outside of the domain fl, one can uniquely determine 
the function u {x,t) for {x,t) £ (M"\r2) x (0,T) . Indeed, this follows from Remark 2.2. Therefore we can 
uniquely determine functions u,9„u at St- Hence, we can replace in p.39p Too with 5*00. Hence, Theorem 
3.1 implies Theorem 3.3. 

Theorem 3.3. Assume that conditions Ii3.34\ l hold. Also, assume that coefficients c(x) ,aa {x) and the 
initial condition fo (x) are so smooth that the solution v of the problem i3.36\) . 113.37^ is such that: 

(a) V & ^ if the function c{x) is unknown, and (h) v £ C^+\'^\ {l^^j if the function Oa {x) 

is unknown. Let the domain fl = {|a;| < R} C M",n > 2. Denote b(x) = l/-\/c (x). Let the function b{x) 
satisfies conditions Ii2.59\) . \2.6U\) . Suppose that conditions of either of Cases 1 or 2 of Theorem 3.1 hold. 
Then conclusions of Theorem 3.1 are true for the inverse problem \3. 32\) - ^3.35\) . 

3.3.2. The second MCIP for the parabolic equation. Two points of Theorem 3.3 arc inconvenient 
ones. First, one needs to reduce the parabolic CIP to the hyperbolic CIP via inverting the transform p. 381) . 
Second, one needs to use a special form of the elliptic operator in p.36p with the restrictive condition (j2.60p 
imposed on the coefficient c{x). Although (j2.60p holds for the case c{x) = 1, still the question remains 
whether it is possible to prove uniqueness of an MCIP for the case of a general parabolic operator of the 
second order. We show in this Section that the latter is possible, provided that one can guarantee the 
existence of the solution of the parabolic PDF for t G {~T, T) and that the function u (a;, 0) is known. 
On the other hand, it was noticed in the paper of Yamamoto and Zou |150j that the measurement of the 
temperature at i = > is often easier to achieve than the measurement at the initial time moment t = 0. 
Thus, that condition likely has a good applied sense. 

Let f2 C M" be either finite or infinite connected domain with the piecewise smooth boundary 9f2, 
F C 057 be a part of this boundary and T = const > 0. Denote F^ = F x {—T,T) . Let L be the elliptic 
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operator in of the form 

Lu= ^ a„ (a;) £>>, (x, t) e Q^, (3.40) 

|q|<2 

a„ e (n) , |a| = 2, aa e C (fl) , \a\ = 0, 1, (3.41) 
MilCI^< XI "«(^)?"^/^2l€I^VCeIR'^Va;e^;//l,^2=cons^. >0. (3.42) 

l"l=2 

The Second Parabolic Coefficient Inverse Problem. Assume that one of coefficients Oq-q (a;) of the 
operator L is unknown in fi, and all other coefficients of L are known in VI. Let the function u G C*'^ ^Qt) 
satisfies the parabolic equation 

ut = Lu + F (x, t) , in . (3.43) 

Determine the coefficient a^g (x) for a; € assuming that the function F (x, t) is known in and that the 
following functions /o (x) ,p (a;, t) and g (a;, t) are known as well 

m(.t,0) = /o(a;),xe 17, (3.44) 

u|p± :=p(a;,<),— lr± ='?(x,i). (3.45) 

Theorem 3.4. Assume that conditions {3.41^ and {S.^i^ are valid and that 

DTfo (x) ^ in n. (3.46) 

Then there exists at most one solution (u, aa„) of the inverse problem (34§-(34§ such that u e C"^'^ (Qt) ■ 
Proof. Assume that there exist two solutions {ui,a^al^ and ^1*2,000''^ Denote u = mi — M2,a = 
Oa} - ai]. Using we obtain 

Ut- L^^'>u = Z{x)D'^''u2, (3.47) 
2(x,0) = 0, (3.48) 

"lr± = |^lr±=0. (3.49) 

Here L^^^ means that the coefficient [x) in the operator L is replaced with a^a} [x) . It follows from p.46p 
that there exists a small number e > such that D°'°U2 7^ in . Denote for brevity g (x, t) ~ D°'°U2 {x, t) . 
By (IMZD 

~ / \ Ut — L^^^u — ± 

a{x) = ,[x,t)& Q 

9 

Differentiating this equality with respect to t and denoting 

h{x,t) = —, (3.50) 
9 

we obtain 

utt - L'^^^ut = h(ut- L^^'>u^ . (3.51) 

We now proceed similarly with the proof of Theorem 3.1. Denote 

V = ut — hu. (3.52) 
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Solving the Ordinary Differential Equation p.52p with the zero initial condition (|3.48p and taking into 
account (|3.50p . wc obtain 



u {x, t) 



K (x, t, t) V {x, t) dr, K {x, t, r) = 7-77-^, {x, t) G 





Next, using (|3^ . (|33T|) and (|3?52|) . wc obtain similarly with (|3TT9|) and ((3?20| 



(3.53) 



r(i) 



< M 



|Vu| (x, t) + \v\ {x, t) + sgn (t) / (|Vw| + \v\) (x, r) dr 
dv 



dn 



.±= 0, 



(3.54) 



(3.55) 



where Lq^'' is the principal part of the operator L^^\ sgn (t) = 1 if t > and sgn (t) = —1 ii t < 0, and the 
positive constant M is independent on x,t,v. Without loss of generality we assume that 



r = {xi = o,|y| <r}. 



(3.56) 



Otherwise we can still obtain (|3.56p for at least a piece of F via a change of variables. Consider such a point 
ccq G f2 that the straight line which is perpendicular to {xi = 0} and passes through xq, intersects F, and the 
segment of this straight line which connects xq and F, lies inside of il. Without loss of generality we assume 
that 



Xq e Ixi > 0,xi + <rj-a 



(3.57) 



where numbers a, rj d (0, 1) , a < 77 were defined in subsection 2.3. 

Let tp{x,t) be the function defined in p.24|) . (|2.25p . In ()2.24p choose the parameter T so large that 
e^/ (2r^) < rj — a. Hence, the domain G,, C Qf, where was defined in subsection 2.3. Square both sides 
of inequality p.54[) . multiply by the function cp'^ {x, t), integrate over the domain G,,. Following Remark 2.2, 

we use the Carleman estimate of Theorem 2.3, assuming that it is valid for the operator L^"^ . In addition, 
we use p.55p . Also, fix the parameter v :— {K, Gr/, /ii, IJ-2) ■ Wc obtain with a new constant M 



C\ J {Vv f ip'^dxdt + CX^ J v^ip^dxdt <M J \^iVvy + w^J ip'-'dxdt + M 

Gjf Gjj Gj^ G 

+CX^ exp [2X7]-"] J (Vf ) 
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(|Vw| + \v\)dT 



ip^dxdt 



dS. 



Using Lemma 3.1, we obtain similarly with the proof of inequality p.23p with a new constant M 

2 



M 



\Vv\ + \v\)dT 



ip^dxdt < 



M 

X 



tp^dxdt. 



Hence, choosing sufficiently large Ai = Ai (Af, K, G^, /i^, fJ-2) > Aq, we obtain 



^ip'^dxdt <CX^exp[2XT]-''] J [(Vu) 

32 G, 



2 , 2 

V 



dsyx > Ai. 



(3.58) 
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Let (5 e (0,1] ~ a) be so small that the point (a;o,0) £ G^-a. Then making p.58p stronger, we obtain 



A exp 



2\{ri-5y 



''dxdt < CX^ exp [2X7]-"] 



dS,\fX > Ai. 



Gn-5 



02 G„ 



Dividing this inequality by A^ exp 2A (rj — 6) 



, we obtain 



J w^dxdt < Cexp|-2A (rj - Sy" - rj^" } ,VA > Ai. 



Gn-s 

Setting here A — cxd, we obtain 



v^dxdt — 0. 



Hence, v{x,t) = in G,i-s- Hence, (|3.53p implies that u{x,t) = in Gn-s- Substituting this in p.47p 
and using D"°U2 {x,t) ^ in Gj^-s, we obtain a (a;) = in G,^-s n {t = 0} . In particular a(xo) = 0. It is 
clear that rotating and moving the coordinate system, one can cover the entire domain VL by paraboloids 
like Gr^-s n {i = 0}. Thus, = 0. This, ([3^491) . Theorems 2.3, 2.4 and Remark 2.2 imply that 

ui {x,t) — U2 {x,t) in Q^. □ 

3.3.3. An MCIP for a parabolic equation with final over determination. Let Dj 

(0,T) and L be the elliptic operator in R", whose coefficients depend only on x, 



Lu = ^ aa (x) D'^u, 

\a\<2 

a„ e C^+f (R"),/3e (0,1), 
Mi|e|'< fla (a^) r <A*2iei',Va;,eeR", Ail, A*2= const. >0. 

\a\=2 



Consider the following Caucliy problem 



ut = Lu in D'^^, 



Jrj, 



(3.59) 

(3.60) 
(3.61) 



(3.62) 
(3.63) 



It is well known that the problem p.62p . p.63p has unique solution |104j . 

The Parabolic Coefficient Inverse Problem with Final Overdetermination. Let To G (0,T) 
and il C R" be a bounded domain. Suppose that the coefficient Uao i^) of the operator L is known inside of 
n and is unknown outside of fl. Assume that the initial condition / (x) is also unknown. Determine both 
the coefficient aag (x) for x G R"\J7 and the initial condition / (x) for x G R" , assuming that the following 
function F (x) is known 



F{x) = u{x,Tq) ,x e 



(3.64) 



Theorem 3.5. Assume that conditions i3.5y\} - ![3.61\} hold, all coefficients of the operator L belong to 
C°° {n), and 



D'^'Tix) ^ 0, in W^\n. 



Then there exists at most one pair vector function (u (x, t) , (x) , / (x)) satisfying conditions i3.6S^) - (3.64^ . 
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Proof. Consider the solution of the foUowing hyperbohc Cauchy problem 



Vtt=Lv in 
v{x,Q)^Q,Vt {x,0)=f{x). 

Then by p.38p u ~ TZv. Hence, for any x G M" the function u {x, t) is analytic as a function of the real 
variable t > 0. We now show that the function u {x, t) can be uniquely determined for (x, t) G x (0, T) . 
Since all coefficients aa € (51), then u £ C°° (51 x (0,T)) , see the book of Friedman [54]- Hence, using 
(|3:62l) and t^M^ . we obtain 

D^+\ {x, To) = [F (x)] , .T e O, fc = 0, 1, . . . . 

This means that one can uniquely determine all t derivatives of the function u (x, t) aX t = Tq for all x G Jl. 
Hence, the analytieity of the function u (x, t) with respect to t implies that this function can be uniquely 
determined for (x,i) e 51 x (0,T). Next, applying Theorem 3.4, we obtain that the coefficient (x) is 
uniquely determined in the domain M"\r2. The classical theorem about the uniqueness of the solution of the 
parabolic equation with the reversed time |541 1112j implies that the initial condition / (x) is also uniquely 
determined. □ 

3.4. An MCIP for an elliptic equation. We now consider an elliptic analog of the Second Parabolic 
Coefficient Inverse Problem. Let 51 C M" be either finite or infinite domain with the piecewise smooth 
boundary 351 and let F C 951 be a part of this boundary. Let T = const > 0. We keep notations of Section 
3.3. Let L be an elliptic operator in 51, 

Lu = ^ aa (x)il>,x e 51, (3.65) 

l«|<2 

a„ e (51) , |a| = 2; a„ e C (U) , |a| = 0, 1, (3.66) 
Mil?l^< ^ a„(x)C"< /l2l^l^V?eR",Vx e51;Ati,Ai2= const. >0. (3.67) 

l"l=2 
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Coefficient Inverse Problem for an Elliptic Equation. Let the function u £ C [Qtj satisfies 
the following conditions 

utt + Lu = F{x,t) in (3.68) 

u(x,0) = fo (x) in 51, (3.69) 



\^±=p{x,t), |^|p± =g(x,t). (3.70) 



Assume that the coefficient a^g (x) of the operator L is unknown in 51 and all other coefficients arc known. 
Determine the coefficient (x) from conditions p.68p - p.70p . 

Theorem 3.6. Assume that fo (x) 7^ in 51 and conditions (3. 65\) - f3.67\ ) hold. Then there exists 

at most one pair of functions (aqq (x) ,u(x,t)) satisfying i3.68\) - (3. 70| ) and such that 

It follows from Theorem 2.4 and Remark 2.2 that the proof of this theorem is completely similar with 
the proof Theorem 3.4. Therefore, we omit this proof here. 

4. Published Results About BK. Given a significant number of publications about BK, it would 
be quite space consuming to tell details about the topic of each one. Therefore, the author provides short 
comments about cited papers. An interesting reader is referred to the corresponding paper for detail. Many 
works cited in this section are devoted to either Lipschitz or Holder global stability estimates for MCIPs. 
Holder stability estimate means ||ai — 02!!^ < C'||i^i — ^2||y ,ct — const. G (0, 1] . Here ai and 02 are two 
unknown coefficients, corresponding to the data Fi and F2 respectively and X and Y are corresponding 
Banah spaces. The case a = 1 is called Lipschitz stability estimate, which is obviously stronger than a < 1. 
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This is why the Lipschitz stability estimate is usuahy much harder to prove than the Holder stability. On 
the other hand, it was briefly mentioned in earlier papers of the author [TSl [55] that usually the Holder 
stability follows almost immediately from BK. To do this, one needs to combine BK either with Theorem 
2.1 or with similar theorems of Chapter 4 of the book |112j . All stability estimates for finding coefficients 
mentioned in this section are conditional stability estimates, as it is usually the case in the theory of Ill-Posed 
problems [Sj I14| I49| 1701 1145) . In other words, some a priori upper bounds are imposed either on certain 
norms of coefficients of interest or on certain norms of solutions of corresponding PDEs. The constant C 
depends on these bounds. In some works the problem of finding an unknown coefficient is replaced with the 
problem of finding the function / {x) in the source term / {x) P {x, t) , where the function P {x, t) is known. 
This problem is almost equivalent to a corresponding MCIP, since the interpretation of / [x) in this case is 
/ {x) = fli [x) — 02 (x), see p.7p and (j3.10p . Since the problem of finding the source term is linear, unlike an 
MCIP, then a priori bounds depend on the function P{x,t). 

4.1. MCIPs for hyperbolic PDEs. Lipschitz stability is established for some of these MCIPs. This 
became possible because the hyperbolic PDE can be solved in both forward and backward directions of 
time. The idea behind Lipschitz stability estimates is to use various combinations of BK with the Lipschitz 
stability estimate for the Cauchy problem for the hyperbolic equation with Dirichlet and Neumann data 
given at the lateral boundary of the time cylinder. For the latter see Theorem 5.1 in Section 5.1 as well as 
originating papers of Klibanov and Malinsky |5T] and Kazemi and Klibanov [71| . 

For the first time, the Lipschitz stability estimate for an MCIP with single measurement data was 
obtained by Puel and Yamamoto |131| . The initial boundary value problem in |131j is 



The MCIP of |131) consists in finding the source function / (x) , assuming that the normal derivative g (x, t) = 
dnU |sj,is known and P (x, 0) ^ in Vl. Hence, this MCIP can be obtained from the MCIP of finding 
the unknown coefficient p {x) by assuming that / (x) = pi (x) — p2 (x) , where pi and p2 are two possible 
coefficients. In this case in (j4.1[) p — p2, u ^ ui — U2, P ~ ui, where ui and U2 are solutions of the problem 
(|4.1[) with p := pi and p := p2 respectively with the same initial conditions 



and the same Dirichlet boundary condition ui \st— ^2 \st— P \st - It was shown in |131] that if f2 = 
{\x\ < R} and T > 2R, then 



Isakov and Yamamoto |64j and Imanuvilov and Yamamoto |601 161] have obtained various Lipschitz stability 
estimates for MCIPs for hyperbolic PDEs with the principal part of hyperbolic operators df --A. Bellassoued 
[23] has proved Lipschitz stability for an MCIP for the equation uu ~ (? ix) Au, where A is a self-adjoint 
elliptic operator of the second order with a;— dependent coefficients and c (x) is the unknown coefficient. 

Imanuvilov and Yamamoto |62| have considered the case of the determination of the coefficient p (x) in 
equation 



A new element here, compared with Theorem 3.1, is that both the function p and its first derivatives are 
involved in equation (|4.2p . On the other hand, the machinery of Theorem 3.1 would require to consider 
derivatives Px^ as "independent" functions. In turn, this would require to use n + 1 independent initial 
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utt ^ Au- p{x)u- f {x) P {x, t) , {x, t) e Qt, 
u (x, 0) = ut (a;, 0) = 0, 
u I = 0. 



(4.1) 



Ui (x, 0) = U2 (x, 0) = P (x, 0) , Uxt {x, 0) = U2t {x, 0) 




Utt = div {p {x) Vu) . 



(4.2) 



conditions. However, only one initial condition was used in |62j . Holder stability estimate was obtained in 
[52] , This result was extended by Klibanov and Yamamoto [55] to the Lipschitz stability. The method of 
[55] is different in some respects from the one of (55], because in [55] a combination of ideas of Theorem 3.1 
and Theorem 5.1 was used. 

Publications cited above in this section used the assumption that the Dirichlet boundary condition 
is known on the entire boundary St and the Neumann boundary condition is known on a part of the 
boundaryc^ri satisfying an appropriate geometrical condition. To obtain logarithmic stability estimates for 
the case when the Neumann boundary condition is known on an arbitrary piece of the boundary, Bellassoued 
[53] has proposed to use the so-called "Fourier-Bros-Iagolnitzer integral transformation" (FBI) with respect 
to t . FBI transforms the hyperbolic equation in the elliptic one, where the operator df — A is replaced with 
df + A. While in [23] this was done for equation utt = Am — p (x) u, Bellassoued and Yamamoto [24] have 
extended this result to the case of equation (|4.2p . We refer to Robbiano |133[ I134j for the introduction of 
the FBI transformation. 

Liu and Triggiani |118| have considered the hyperbolic equation with the damping term q (x) ut. 



The MCIP of |118j consists in determining the coefficient q (x) from the Dirichlet data u \rj,= f {x, t) , where 
r C 9f2 is a part of the boundary 9f2 satisfying an appropriate geometrical condition. The Lipschitz stability 
estimate for this MCIP was obtained in |118] . 

As to other publications about the use of BK for MCIPs for hyperbolic PDEs, we refer to works of 
Khaidarov [73] , Isakov [65] IMl [62] , Doubova and Osses [47] , Baudouin, Crepeau and Velein [10] , Baudouin, 
Mercado and Osses [8] , Yuan and Yamamoto |152j and Liu and Triggiani [119, 120| . All above results for 
MCIPs for hyperbolic PDEs are obtained under assumptions like the one in (|2.60|) , which is imposed on the 
coefficient c {x) in the principal part of the hyperbolic operator. This is because conditions like (|2.60p are 
the only known ones which guarantee the existence of the Carleman estimate for the hyperbolic case on the 
one hand, and can be directly analytically verified for generic functions on the other hand. Note that (|2.60|) 
is valid of course for the case c{x) = 1. 

4.2. MCIPs for parabolic PDEs. A comprehensive survey about the topic of this section as well as 
about some related topics can be found in the paper of Yamamoto |151j . Thus, the author refers to |151| for 
further references. It is a long standing well known open problem to prove uniqueness theorems for MCIPs 
for parabolic PDEs with single measurement in the case when the regular initial condition is given at {t = 0} 
and the equation is valid only for t € (0,T) . This is why the only case when the uniqueness can be currently 
proven for this type of data is the one of Section 3.3.1, where the inverse Reznickaya's transform ()3.38p was 
used to obtain the MCIP for a similar hyperbolic PDE. Hence, we discuss in this section only the case when 
the equation is valid for t S (— T, T) and the data are given at {t = 0} as well as on at least a part of the 
lateral boundary. The only exception is the case of nonlinear parabolic PDEs. 

An inconvenience of the conventional Carleman estimate of Theorem 2.4 (Section 2.3) is that it is valid 
only in the paraboloid G^, which is a subdomain of the time cylinder Q^. On the other hand, Fursikov 
and Imanuvilov |551 158j have proved a radically new Carleman estimate for an arbitrary parabolic operator 
of the second order. This estimate is valid in the entire time cylinder Q^, although the Carleman Weight 
Function exponentially decays to zero at t — >■ ±T. Using this fact, Imanuvilov and Yamamoto |59j have 
proved, for the first time, the Lipschitz stability estimate for an MCIP for a general parabolic equation with 
cc— dependent coefficients. In [59j the Dirichlet boundary condition for the solution of the forward problem 
was known at the entire boundary, whereas the Neumann boundary condition was known on any piece of 
the boundary. Starting from [59j , Carleman estimates of the Fursikov-Imanuvilov type became popular in 
the inverse problems community, see, e.g. the papers of Baudouin and Pucl [7] and Cristofol, Gaitan and 
Ramoul [IS]. 



iLtt = Au + q (x) ut, (x, t) £ Qt, 



u l^x, —j = uo {x) , ut l^x, — 

dnU I St = 5 {x, t) . 
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Yamamoto |151) has obtained the Holder stabihty estimate for the case of the equation Uf ~ div (p (x) Vw) 
with the unknown coefficient p (x) in the case when two boundary conditions are given at any part F of the 
boundary 917. The main difference between this resuU and the one of Theorem 3.4 is that the machinery 
of Theorem 3.4 would require to treat first derivatives of the function p (x) as independent functions, which 
would lead, in turn to the necessity to use {n + 1) conditions at {t ^ 0} . Unlike the latter, only one condition 
at {t = 0} is used in |151j . It was pointed out in Remark on page 41 of |151j that if one boundary condition 
would be known at the entire boundary dCl and the second one would be known only at F, then the Lipschitz 
stability estimate would be obtained. In the latter case the Carleman estimate of the Fursikov-Imanuvilov 
type would be used. 

Yamamoto and Zou |150j have considered an MCIP for the equation 



Let uj d Q he subdomain of the domain fJ. An interesting new feature of this work is that the inverse 
problem consists in the simultaneous reconstruction of both the coefficient p (x) and the initial condition 
fi (x) = u (x, 0) , assuming that the following functions / and g arc given 



where — const, e (0,T). First, using the technique of [59], they proved Lipschitz stability estimate for 
the function p (x) . Next, they proved logarithmic stability for the initial condition ^ (x) using the method 
of logarithmic convexity of Payne |129| . This method works for the case when the corresponding elliptic 
operator is self-adjoint. In addition, they have constructed a numerical method, which is based on the 
minimization of the Tikhonov functional. A careful convergence analysis was provided. That analysis was 
confirmed by a number of numerical experiments. 

Surprisingly, Lii |122| . has applied BK, for the first time, to an MCIP for the stochastic parabolic PDF. 

The assumption of Theorem 3.4 that coefficients of the operator L are independent on t was imposed 
only for the sake of simplicity. In fact, one can allow all coefficients, except of a^o [x) , to be dependent 
on both X and t. A direct analog of Theorem 3.4 is valid in this case. To prove it, one should use the 
Carleman estimate of Theorem 2.5 (for ig^''), a certain change of variables and the assumption that the 
target coefficient aa^ (x) is known for x G F, see works of the author [75j [77l |82] as well as Theorem 1.10.7 
in the book [T3]. The same is true for Theorem 3.6. This idea was used in the works of the author discussed 
in the next paragraph. 

BK was also applied to MCIPs for nonlinear parabolic PDEs. In the Id case the author has considered 
the inverse problem for the equation 



where du^^F (y) ^ 0,Vy € ]R"+'^ [T^- Let {xi}"^-^^ C (0, 1) be a sequence of points, such that Xi ^ xj \ii ^ j. 
The inverse problem in |79| consists in determining the vector function q (u) ~ {qi, qn) (u) , assuming that 
the functions fi (t) — u {xi,t) ,i £ [l,n + 1] are known. It was assumed that Ux > in. [0, 1] x [0, T] . This 
inequality can often be established via the maximum principle. Uniqueness theorem was proved in [79j . The 
first step was to introduce a new spatial variable z and a new function v (z, t) via u (v (2, t) , t) := z. 
In [85] and Chapter 4 of [86] the author considered an MCIP for the nonlinear parabolic equation 

Ut=F {Ux^a:i,-,Ux^x„,Ua:^, u^^ , u, x,t, q {u, X2, ...,x„)) , X G {xi € (0, 1) , y G O'} , t G (0, T) , 

where y = (x2, ...,x„) , and du^.^.F{z) G [^1,^2] i^-^^ € R^"+^, where fJ-i, fJ-2 = const. > 0. Here fl' C R"~i 
is an arbitrary domain. The following functions ipi, ip2,ipi, 1P2 were given in [85| 




<Pi (y, t)^u (0, y, t) , -01 (y, t) = (0, y, t) , ip^ (y, t) = u (1, y, t) , -02 (y, t) = Uxi (1, y, t) , (y, t) e n' x (0, T) . 



lit = Alt + p (x) u, (x, t) G Qt, 

U I St =">! i^it) 



u {x, e) = f {x) ,u |ac^x(o,T)= g {x, t) , 



ut= F {u^x, u^, u, qi (u) , 



qn{u)),{x,t) G (0,1) X (0,T), 
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It was required to reconstruct the function g (u, 2:2, a;„) . The first step was again to introduce a new 
spatial variable z and a new function v{z,y^t) via u{v (z^y,t) ,t) := z. Next, uniqueness theorem was 
proved. However, a stability estimate was not established in |85j . Furthermore, unlike the linear case. Holder 
stability estimate does not follow automatically from BK in this nonlinear case. Thus, the Holder stability 
estimate for a similar inverse problem was proved in the paper of Egger, Engl and Klibanov j48j for the case 
of the equation 

■Ut =Uxx + q [u) + / (x, t) 

with the unknown function q. In addition, a numerical reconstruction procedure was developed in |48j 
via minimizing the Tikhonov functional. Numerical results were also presented in [48]. Other MCIPs for 
nonlinear parabolic PDEs were treated via BK in Boulakia, Grandmont and Osses [21] and Kaltenbacher 
and Klibanov [71]. 

Some other MCIPs for parabolic PDEs were treated by various modifications of BK in papers of Bel- 
lassoued and Yamamoto [25l, Benabdallah, Dermenjian and Le Rousseau ^H], Benabdallah, Gaitan and Le 
Rousseau [29] and Poisson (130| . 

Isakov [SS1[S7] has proved uniqueness of the parabolic MCIP with the final overdetermination without the 
assumption of Theorem 3.5 of the knowledge of the target coefficient in the domain Vt. The initial condition 
is also known in [65] [67] . BK was not used in [65] [67] . 

4.3. MCIPs for the non-stationary Schrodinger equation. Baudouin and Puel [7] were the first 
ones who has applied BK to the MCIP for the non-stationary Schrodinger equation 

iut + Au + 5 (x) u = 0, (x, t) e Qt, i"^ = —1, 

u{x,Q) = uo{x), (4.3) 

U I St =h{x,t) . 

The inverse problem in [7] consists in determining the coefficient q [x) from the Neumann boundary condition 

5„u|r^=.9(.T,i),rT = rx (0,T), (4.4) 

where P C dU, is a part of the boundary satisfying an appropriate geometrical condition. First, following 
the idea of [59], an analog of the Carleman estimate of Fursikov and Imanuvilov [55] [58] was proved. Next, 
the Lipschitz stability for the MCIP ([4?3| . ([Q]) was established. 

For three other results of the topic of this section we refer to works of Baudouin and Mercado [9], 
Mercado, Osses and Rosier |123j and Yuan and Yamamoto |153j . 

4.4. Non-standard PDEs. Baudouin, Cerpa, Crepeau and Mercado [TT] have considered the Kuramoto- 
Sivashinsky equation (KS) in Id with Qt = (0, 1) x (0,T), ^ 

Ut + (cr {x) Uxx)xx + 7 (a;) Uxx + u-Ux= g, (4.5) 

u{x,0) =uo{x), (4.6) 



u\x=o = hi {t) , u \x=i^ h2 (t) , (4.7) 
Ux I 1^=0 = /^3 (t) , Ux U=i= /14 (t) ■ (4.8) 

Note that KS is a nonlinear equation. Conditions (|4.5|) - ()4.8|) define the forward problem for KS. First, an 
existence, uniqueness and stability theorem for this problem was proved in [11]. In the inverse problem the 
following functions were assumed to be known 

u {x, To) , Uxx (0, t) 

(0,0, 

where To ~ const. G (0, T) . As to the knowledge of the function u (x, Tq) , see the arguments in the beginning 
of Section 4.2. The Lipschitz stability for the inverse problem was proved in [TT] . 
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Cavaterra. Lorenzi and Yamamoto considered an MCIP for a PDE whose principal part was a 
hyperbolic operator — [x) A and lower order terms included Voltcrra integrals 



t 




(4.9) 







The Lipschitz stability estimate was obtained in [43]. Romanov and Yamamoto |139j obtained the Holder 
stability estimate for an MCIP for a hyperbolic- like PDE with integrals (|4.9|) . We also refer to the work 
of Buhan and Osses |53] , where logarithmic stability estimate for a hyperbolic- like coupler system of PDEs 
with integrals like the one in (|4.9p was obtained, a numerical method, based on the minimization of the 
Tikhonov functional, was developed, and numerical results were presented. 

4.5. Coupled systems of PDEs. In this section we do not present coupled systems under discussion 
because they are well known. Still, we present two systems of PDEs, which are not conventionally known. 
An MCIP for the Maxwell equations was considered by the author in [78| . Let e (x) and a {x) , a; G be 
the dielectric permittivity and the electric conductivity coefficients respectively. It was required to find both 
of them simultaneously, given the magnetic vector field 'H.{x,t) outside of the domain of interest. It was 
assumed that the magnetic permeability coefficient ^ {x) = 1. Uniqueness theorem was proved. We also refer 
to the book of Romanov and Kabanikhin fl36j for inverse problems for the Maxwell's system with impulsive 
sources. 

Yamamoto [149] studied an inverse source problem for the Maxwell's system. In this case the unknown 
source vector function f (x) is three dimensional. Both electric and magnetic vector fields were known at the 
boundary. Uniqueness theorem was proved. We also refer to the papers of Li |113j and Li and Yamamoto 
|114| for the cases of bi-isotropic and anisotropic Maxwell's system respectively. In |113j the Lipschitz 
stability estimate was proved, and in the paper |114| the Holder stability was establishes. Note that in |114j 
the unknown coefficients are actually matrices of dielectric and magnetic permeability coefficients which are 
independent on one spatial variable but dependent on time t. 

Bellassoued, Cristofol and Soccorsi [27] obtained Holder stability estimate for an MCIP for the Maxwell's 
system. In this case both functions e (x) and fi (x) were unknown and a (x) = 0. Tangential components 
of both magnetic and electric field were measured at the boundary for two sets of initial conditions, i.e. 
measurements were conducted twice. 

Imanuvilov, Isakov and Yamamoto [63] considered the MCIP of the reconstruction of all three elastic 
coefficients A {x) , fi (x) , p (x) in the time dependent Lame system. Two sets of initial data were used and 
they have generated two sets of boundary conditions, which were used as the data for MCIP. The Holder 
stability estimate was obtained. 

In |63j the data for the inverse problem were given at the entire boundary. Unlike this, Bellassoued, 
Imanuvilov, and Yamamoto [26] considered the MCIP of recovering of elastic coefficients A (x) , // (x) , p (x) 
of the time dependent Lame system in the case when the data for the MCIP are given at an arbitrary part 
r C dn of the boundary for t G (0, T) . In other words, the previous idea of Bellassoued [23] and Bellassoued 
and Yamamoto [21] (Section 4.1) was extended from a single hyperbolic equation to the to the case of 
Lame system. Similarly with [331 [21] the Fourier-Bros-Iagolnitzer transformation was applied. Logarithmic 
stability estimate was obtained in [26]. 

Liu and Triggiani |117| considered the following 2x2 coupled system of Schrodingcr equations in Qt 



iwt + Aw 
izt + Az 



m {x) ■ Vw + n{x)w + q {x) ■ Vz + p {x) z, 
{x) ■ + a {x) z + ij] {x) ■ "S/w -|- q {x) w, 



( ^ 




wo {x) ,z (x, ^ ) = zo {x) , 




dnW I St = 51 (a;) ,9„Z |sr= 52 (2;) • 



The MCIP in |117j consists in the determination of the pair of unknown coefficients {p, q) {x) ,x £ from 
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the Dirichlet boundary data 



w |rT= fi , z Ir^^ /2 {x,t) , 

where T C dQ is a part of the boundary satisfying an appropriate geometrical condition. Uniqueness theorem 
was proved in I117j . 

Fan, Di Cristo, Jiang and Nakamura [50] proved Lipschitz stabihty estimate for an MCIP for Navicr- 
Stokes equations. The MCIP consists in recovering of the viscosity function from a single boundary mea- 
surement. 

Wu and Liu [148] considered the MCIP for the thcrmoelastic system with memory (also, see their 
preceding work |147| V That system is 

Utt — aAu — div U+7VW = a {x, t) p (x) , (a;, t) e Qt, 

t 

— J k{t — T) Av (x, T)dT + ^ div Ut = 0, (x, t) e Qt, 


u {x, 0) ~ Uq (x) , Ut {x, 0) ~ Ui (x) , V (x, 0) — Vi (x) , 
U I St = |St= 0. 

Here u = (ui, U2, 1*3)"^ and v are displacement and temperature respectively. Let w C fi be a subdomain. The 
MCIP in 1 148] consists in the recovery of the source term p (x), given the vector function f (x, t) = u (x, t) 
for (x, t) € Lo X (0, T) . The function a (x, t) is assumed to be known. Lipschitz stability estimate for this 
MCIP was obtained in [T48] . 

Cristofol, Gaitan and Ramoul [ISj considered the following parabolic 2x2 system 

Ut = Au + a{x)u + b (x) v, (x, t) G Qt, 
Vt = Av + c (x) u + d (x) w, (x, t) e Qt, 

U (x, 0) = Uq (x) , V (x, 0) = Vq (x) . 

This is a reaction-diffusion system. Let w C fJ be a subdomain, to E (0,T) and T' = (fo + T) /2. The data 
for the MCIP in [55] are the following functions 

Au{x,T') ,u{x,T') ,v{x,T') ,vt Lx(to,T) ■ 

The MCIP consists in determining the vector function (6,ito,wo) .The function b can be replaced with the 
function a. A new point of [45j is that only one function Vf is measured in cj, unlike the conventional way 
of measuring both functions u and v. The Lipschitz stability estimate for the coefficient b (x) as well as the 
logarithmic stability estimate for initial conditions uq (x) ,vo (x) were obtained in [45] . 

5. Stability Estimates for Hyperbolic Equations and Inequalities With Lateral Cauchy 
Data and Thermoacoustic Tomography. It was shown in Section 2.2 how Carleman estimates lead to 
Holder stability estimates for ill-posed Cauchy problems for PDEs, including inequalities. In this section 
we obtain Lipschitz and logarithmic stability estimates for the case of hyperbolic PDEs. We also show how 
these estimates help to specify QRM for the hyperbolic case. In Section 5.5 we cite related published results. 

5.1. Lipschitz stability. As usual, we consider the simplest case when the domain of interest 17 is a 
ball, fl = {\x\ < R} C M", where R — const. > 0. Let functions c (x) , bj (x, t) satisfy the following conditions 

c (x) G (H) , c (x) 7^ in n, (5.1) 

5,GC(Q^),j=0,...,n+l. (5.2) 
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Let the function u £ C [Qxj ^-nd the function f & C [Qxj ■ Suppose that the function u is a solution of 
the following hyperbolic equation 

n+l 

utt = <? (x) Am + ^ bj {x, t) + 60 (a;, t)u + f (x, t) , {x, t) e (5.3) 

i=i 

where Un+i Ut- Consider Dirichlet and Neumann boundary conditions for the function u at the lateral 
side of the time cylinder , 

u\g±= p{x,t) ,dnu\g±= q{x,t) . (5.4) 

Problem 5.1. Given conditions (|5.3p . (|5.4p . estimate the function u £ C"^ {Q^^ in the time cylinder 
via functions p, q and /. 

The of this section works for a more general case of a hyperbolic inequality. Specifically, we consider the 
following problem. 

Problem 5.2. Let A = const. > 0. Let the function m e (Qt^ satisfies the following pointwise 
inequality in the cylinder 

\uu-c^x)Au\ < Ai\Vu\ + \ut\ + \u\ + \ f\) y ix,t) e Q^, (5.5) 

where the function / e L2 {Qx) ■ Estimate the function u via functions p, q and /. 

Since Problem 5.2 is more general than Problem 5.1, we study only Problem 5.2 in this section. Theorem 
4.1 provides the Lipschitz stability estimate for Problem 5.2. We now reformulate condition (|2.60p for the 
case of the operator df — {x) A in a stronger form as 

(x, Vc"^ (x)) > a = const. > in H, (5.6) 

where a > is a certain number and (, ) is the scalar product in R". Hence, there exists a sufficiently small 
number £ — e (^a, || G (0,i?) such that 

(x- xo, Vc~^ (x)) > ^ > in n^xo e {|a;o| < e}. (5.7) 

Inequality (|5.7p guarantees the Carleman estimate for the operator df — c^ {x) A (Theorem 2.5). This estimate 
is also guaranteed if c = const. ^ 0. 

Theorem 5.1. Let the domain ~ {\x\ < R} . Let a > and d > I be certain numbers. Let 

c e C^n) ,c-^x) e [l,d],||Vc-2||^^^_^ <d. (5.8) 

Ln the case c 7^ const, we assume that condition i5.6] } is fulfilled. Let the function u G {Qt^ satisfies 

inequality I15.5\) with Dirichlet and Neumann boundary conditions ( [5.^[ ). Then there exists a constant 7/q = 
rjQ (i?, d, a) G (0, 1] such that if 

T>^, (5.9) 



then with a constant K = K [A, R, T, d, a) = const. > the following Lipschitz .stability estimate holds for 
the function u, 
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In particular, if c(x) = 1, then one ean take rj^ = 1 and in i5. 9\) T > R. 

Proof. In this proof K ^ K (A, R, T, d, a) denotes different positive constants depending on listed 
parameters. We note first that the constant P {xq, fl) in (|2.6ip can be estimated as P {xq, fi) < 2R. Hence, 

we can set in Theorem 2.5 ?7o = Vo (^Ri \\'^^'^\\ci(p'^^ ^ C-*' ■ '^^ can choose a sufficiently small 

e in (|5.7p and then choose 77 such that 



{R + s 



,2 



r y2 ,^0 c(0,l), (5.II) 



Similarly with Section 2.4 let 7 = const £ (0,£^/9) , 



C (x, t) = \x\^ - j]t^, if {x, t) = exp [A^ {x, t)] , (x, t) e 



T' 

= (x, i) > 7, X e f7} = I |xf - 7]t^ > 7, .T e 0} , 



where A > 1 is a large parameter which we define later. By (|5.1ip 

G^c{\t\<T}. (5.12) 
Choose a sufficiently small number S such that 7 + 3(5 G (O, £^/9) . Hence, 

^7+35 7^ and C G'-y_|-25 C C G-,,. (5.13) 

Introduce a function G C*^ (Qt^ satisfying 



1, (x, t) e 6*7+25, 

X5 (x, i) = <( 0, (x, t) e Q±\G7+5, (5.14) 
between and 1 otherwise. 

The existence of such functions is well known from the Real Analysis course. Let 

v{x,t) =u{x,t)xs{x,t) . (5.15) 

Multiplying both sides of ([53]) by Xs and using (fO| . ()5.12p . (|5.14p and (|5.15p . we obtain 

- c2 (x) Ai;| < If (|Vt.| + + + |/|) + if (1 - Xs) (|Vii| + + \u\) , (x, t) e G^, (5.16) 

ls±= Xsf^dnv |_g±= Xi''/' + <^9„X5- (5.17) 
Squaring both sides of (|5.16p and using Theorem 2.5, we obtain 



K 



> GA (^iVwr + vtj + GA^wV + div t/ + Vt, in G^, VA > Ao. 

where the vector fimction {U, V) satisfies conditions (|2.62p . (|2.63p with the replacement of u by v. Hence, 
(|5.14p implies that U ^ V ^ on {{x,t) : ^ {x,t) ^ j, x £ fl} . Hence, integrating the letter inequality over 
G-y and using Gauss-Ostrogradsky formula and (|5.17p . we obtain 

G y A (\Vvf + v^^ if^dxdt + GA^ J v^ip^dxdt <K J (jVvf + + + g^^ ip^ 
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Let Ao > 1 be the number of Theorem 2.5. There exists a number Ai = Ai(Ao,C, A') > Aq such that 
CA/2 > A',VA > Ai. Hence, 



dxdt < 



< Kexp[2A(7 + 2(5)] \\u\\lur,±\ + Ke^^"' 



lbllHi(s±) + ll9llL(s±) + ll.9llL(s±) 



By ([CT^ and (jETi)) 



\ j {\Vv\^ +vl +v^^(p^dxdt>\ j {\Vv\^ +vl +v^^ (p^dxdt 



Hence, 



\\u\ 



< Kexp {-2XS) ||u||^i(Q±) + Ke^^^ 



>exp[2A(7 + 3<5)]||u||^,(c^^^^). 

lbllffi(s±) + ll9llL(s^) + ll/llL(Q^) 



Note that 



G 



7+35 



n{t = 0} = {x:\x\e (v/7 + 3(5,i?)} D {x : \x\ e 



(5.18) 



(5.19) 



Choose now a pomt such that \xq\ = 3-\/7 + 35. Then |a;o| < £. Consider now an arbitrary point 
ye{\x\< VlTis] . Then 



\y-XQ\> \xa\ - \y\ = 3/7+35 - \y\ > 3^7 + 3(5 - ^7 + 35 = 2/7+35 > ^7 + 35. 



Hence, 



{kl < V7 + 35} c {\y~xo\ > /7 + 35} 



(5.20) 



Introduce now the domain G-y (xq) as 



G. 



(.To) = ||x - xol^ - Tjt'^ > 7, x e ilj 



It follows from (|5.19p and (j5.20p that there exists a sufficiently small number a — a (e) such that 

{t e (0, a)} C U G^+3s (xo)] ■ (5.21) 

Since 7 + 35 £ (0,e^/9) and Ixol = 3^7 + 35, then ((5+T1) implies that G^ {xq) C {|t| < T} . Next, we use 
(|5.7p and Theorem 2.5 and obtain similarly with (|5.18l) 

MIhg^^M) < ^exp(-2A5) + Ke'^'"' [lbll?^.(5±) + Il9lll(5±) + I1/I1L(qJ) 

Combining this with (jS.lSp . we obtain 



< Kexp{^2\S)\\u\\%,^Q±-^+Ke^^^^ 

35 



iipiihi(5±) + ii'^iiL(s±) + ii/iiL(q^) 



This, (j5.2ip and the mean value theorem imply that there exists a number tg G [OjCr] such that 

ll"(2^7^o)llHi(f2) + ll"t (^'^o)lli,(n) 



(5.22) 



< Kexp{-2\S) !|u||^i(Q±) + A'e^^-^ 
Let utt - (a;) u := Z (a;, t) . Then by (|53|) 

|Z| <^(|Vii| + |7.,| + |u| + |/|),V(x,i)e 

Consider the initial boundary value problem with reversed time 

Utt - (? {x)u^Z(x,t) ,y{x,t) <E {x en,te (-T,to)}: 
u (x, to) = Mo (x) , Ut (x, to) = wi (x) , 
u I (a;,t)eaox(-T,io) =p(a;,0 ■ 



(5.23) 



(5.24) 
(5.25) 
(5.26) 



Next, consider the same initial boundary value problem but in the time cylinder (x, t) e {x e f7, < G {to, T)} . 
Recall that the hyperbolic equation can be solved in both positive and negative directions of time. Hence, 
the standard method of energy estimates being applied to two latter problems combined with inequalities 
and leads to 



\\u\\m{Q±) < -f^cxp {^2X6) ||u||^i(Q±) + Ke ' y|p||^i(5±) + 1191^2(5*) + ll/lli2(Q±) 



(5.27) 



Choosing A = X(K,6) so large that Kcxp{—2XS) < 1/2, we obtain the target estimate (|5.10p from (|5.27p . 
□ 

To apply QRM (Section 5.4.2), we need Theorem 5.2. The proof of this theorem is almost identical with 
the proof of Theorem 5.1. 

Theorem 5.2. Let the function u Q (Qt^ satisfies the Dirichlet and Neumann boundary conditions 
^5-4^ as well as the following integral inequality 




ri+l 



Utt 



{x) Au — bj {x, t) Uxj — bo {x, t) u dxdt < S"^ 



and the rest of conditions of Theorem 5.1 is in place. Here S ~ const. > 0. Then there exists a constant 
Va ~ Vo {Ri a) e (0, 1] such that if T > R/ y/rf^, then with a constant K = K {A, R, T, d, a) ~ const. > 
the following Lipschitz stability estimate holds for the function u, 



l'^ilffi(Q±) ^ B 



In particular, if c{x) = 1, then one can take tJq = 1 and T > R. 

5.2. Thermoacoustic tomography. In thermoacoustic tomography (TAT) a short radio frequency 
pulse is sent in a biological tissue, see papers of Agranovsky and Kuchment [5] and Finch and Rakesh [53]. 
Some energy is absorbed. It is well known that malignant legions absorb much more energy than healthy 
ones. Then the tissue expands and radiates a pressure wave. The propagation of this wave can be modeled 
as the solution of the following Cauchy problem 



Utt = (x) Au,x G M^, 
uix,0) = fix),ut (x,0) =0. 



(5.28) 
(5.29) 



The function u{x,t) is measured by transducers at certain locations either at the boundary of the medium of 
interest or outside of this medium. The function / (x) characterizes the absorption of the medium. Hence, if 
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one would know the function / (x), then one would know locations of malignant spots. The inverse problem 
consists in determining the initial condition / (x) using those measurements. Hence, we obtain the following 
inverse problem 

Inverse Problem 5.1 (IP5.1). Let C M"^ be a bounded domain with d^l G C^. Consider the Cauchy 
problem ((5?28)) . (|09)) . Suppose that 

f (x) ^0,x (5.30) 

Determine the function / (x) inside of the domain il assuming that the following function p (a;, t) is known 

u\sT=pix,t). (5.31) 

We show below in this section how Theorem 5.1 implies the Lipschitz stability estimate for IP5.1. 
Assume that c{x) = 1 for x G M^\f2. Then (|5.28p - (|5.3ip imply that the function u{x,t) is the solution of 
the following initial boundary value problem in (M'^Xfi) x (0, T) 

Utt =Au,xe R^\n, t G (0, T) , 
u{x,0) =utix,0) =0,x eR^\n, (5.32) 
U I St =P{x,t) . 

Under certain well known conditions the solution u G ((]R'^\i7) x (0, T)) of this problem exists, is unique 
and depends continuously on the boundary data p {x, t). Just as above consider the case = {|a;| < R} . For 
|x| > i? let p (x, t) ^ p (x, t) \x\ R-^p (x) , where the function p G (|a;| > i?) , p (x) = for x G {|a;| > 3i?} 
and p{x) — 1 for x G G [R, 2R]} . Let v (x, t) = u {x, t) —p {x, t) . Substituting v in equations (|5.32p and 
using the standard method of energy estimates (see, e.g. the book of Ladyzhcnskaya |106| for this method), 
we obtain the following stability estimate |ju||^2((R3\f2)x(o t)) — ^' \\p\\h^{St) ■ ^'^^ ^ {^^t) = (^nU \st ■ Hence, 
trace theorem implies that 

h\\L2{ST) = II^""IIl2(5t) ^ ^ WpWhhSt) ■ (5-33) 

Here C > denotes different positive constants depending only on R and the function p (x) . 

Consider now the even extension u (x, t) with respect to t of the function u (x, t) . Assuming that u G 

(^Qrp) , wc obtain u G {Qt^ a-nd also 

Utt = (x) Au, (x, t) e Q^, 

where functions p and q are even extensions of functions p and q respectively. In addition, by (j5.29p 

u(x,0) = /(x). (5.36) 

Hence, Theorem 5.1, ()5.33p - ()5.36p and trace theorem imply Theorem 5.3. We need here u G H'^ (R^ x (0, T)) 
in order to make sure that the function p G (St) ■ 

Theorem 5.3. Let the domain D, ~ {|x| < R} . Let the function c £ (M.^) ,c = 1 m M^\17 and 
also c satisfies conditions 115. 8\) . In the case c 7^ const, we assume that condition 115. 6\) is fulfilled. Suppose 
that functions c (x) and f (x) are such that there exists the solution u G C"^ (R^ x [0,T]) n iJ"^ (R^ x (0,r)) 
of the Cauchy problem h5. 28\) . i5.29\) . Also, let condition i5. 30\) be satisfied. Then there exists a constant 
Vo ~ '?o(^i'^i'^) S (Oi 1] such that if T > R/ y/ff^, then with a constant K = K {A, R,T, d^ a) > the 
following Lipschitz stability estimate holds for IP5. 1 

\\f\\L.in)<K\\p\\uHST)- 

In particular, if c (x) = 1, then one can take rj^ = \ and T > R. 
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(5.34) 
(5.35) 



5.3. Logarithmic stability in the case of a general hyperbolic operator of the second order 
with a;— dependent coefficients. Condition (j5.6|) is used in Theorem 5.1 because it is linked with the 
existence of the Carleman estimate for the hyperbolic ease, see the end of Section 4.1. Clearly (j5.6p is a 
restrictive condition. Therefore, the next question is whether a stability estimate can be obtained for an 
analog of IP5.1 in the case of an arbitrary hyperbolic operator of the second order. For the first time, this 
question was positively addressed by the author in |95j . We formulate main results of j95j in this section 
without proofs. 

The idea is to apply an analog of the above Reznickaya's transform p.38p . This way the hyperbolic 
PDE is transformed in a similar parabolic PDE. And the function / (x) becomes the initial condition for that 
parabolic PDE. On the other hand, logarithmic stability estimates for the inverse problem of determination 
of the initial condition of a general parabolic equation from lateral Cauchy data were obtained by Klibanov 
[55] in the case of a finite domain 51 C M" and by Klibanov and Tikhonravov [?] in the case of an infinite 
domain C R". Thus, modifications of these results can be applied. Results of both publications [SHI ?] 
were obtained via Carleman estimates. The difference between logarithmic stability estimates for initial 
conditions in [89l ?] and those in the book of Payne [129] is that in [891 ?] the case of a general elliptic 
operator with coefficients depending on (x, t) was considered. On the other hand, the technique of [129] 
works only for self-adjoint elliptic operators with coefficients depending only on x. Since the method of 
[89l ?] works not only for parabolic PDEs but for integral inequalities as well (see Theorem 2.2 in |95] as 
well as Theorem 2.3 in Section 2.2 above), then it enables us to prove convergence of the QRM, unlike the 
technique of [129] . 

We refer to Li, Yamamoto and Zou [115] for another logarithmic stability estimate of the initial condition 
of a parabolic equation with the self-adjoint operator L in a finite domain. A Carleman estimate was also 
used in this reference. An interesting feature of [115] is that observations are performed on an internal 
subdomain for times t € {t,T) where r > 0. In addition, a numerical method was developed in [115] . 

5.3.1. Statements of inverse problems. Let O C {xi > 0} be a bounded domain with the boundary 
dn e C^. Denote P ^ {xi = 0} , Pt ^ P x (0,T) ,Vr > O.Let /c > be an integer and a € (0, 1). Consider 
the elliptic operator L of the second order, 

n n 

Lu = ^ flij- {x) u^^x, + ^ bj {x) + bo (x) u,x e M", 
a,j e C'^+" (R") n (R") , bj,ba G C'=+" (M") , fc > 2, a e (0, 1) , 

n 

Ml \vf < X! - l'7l^ ,Va;,?7 € R"; fj.^^, fi^ = const. > 0. 

Let the function / (x) be such that 

f eCP (R") , p > 3, / (x) = 0, a; e R"\n. (5.40) 
Consider the following Cauchy problem 

utt = Lu,x eW,t e {0,oo) , (5.41) 
u(x,0) = /(x),Ut(.T,0) = 0. (5.42) 

We use everywhere below in Section 5 the following assumption. 

Assumption 5.1. We assume that in (|5.38p . (|5.40p integers k > 2,p > A, coefficients of the operator 
L and the initial condition / are such that there exists unique solution u £ C* (R" x [0,T]) ,VT > of the 
problem (|5.4ip . (|5.42|) satisfying 

ll"llc^(R"x[o,T]) <^e'^^:VT>0, (5.43) 

where the constants B ~ B {L^B) > 0,d ^ d(^L,^ >0 depend only from the coefficients of the operator L 
and the upper estimate B of the norm H/Hj^p^'q) ^ B. 
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(5.37) 

(5.38) 
(5.39) 



Note that (|5.40p as well as the finite speed of propagation of the solution of problem (j5.4ip . (|5.42p guar- 
antee that the function u (x, t) has a finite support ^I* (T) C M", Vi £ (0, T) , VT > 0, see, e.g. the book of La- 

dyzhenskaya [l06]. Hence, (R" x [0,T]) in Assumption 5.1 is actually the space (^~(T) x [0,T]) . Us- 
ing the classical tool of energy estimates |106| , one can easily find non-restrictive sufficient conditions imposed 
on coefficients of the operator L and the function / guaranteeing the smoothness u G (R" x [0, T]) , VT > 
as well as ()5.43|) . We are not doing this here for brevity. We consider the following two Inverse Problems. 
Inverse Problem 5.2 (IP5.2). Suppose that conditions Ii5.37\ )- f5.40\ ) and Assumption 5.1 hold. Let 
u G (M" X [0,r]) ,VT > 6e the solution of the problem \5.41% ^5.4'^ . Assume that the function 
f (x) is unknown. Determine this function, assuming that the following function (pi (x, t) is known 

u\s^='Pi{x,t). (5.44) 

Inverse Problem 5.3 (IP5.3). Suppose that conditions |5.g7[ )-( [5.^0[ ) and Assumption 5.1 hold. Let 

u G (M" X [0,T]) ,VT > Q he the solution of the problem (^4^ . Assume that the function f (x) 

is unknown. Determine this function, assuming that the following function ix,t) is known 

u \xePcs,= V2i^it) ■ (5.45) 

IP5.2 has complete data collection, since the function ip^ is known at the entire boundary of the domain 
of interest ft. On the other hand, IP5.3 is a special case of incomplete data collection, since 51 C {xi > 0} . 

In stability estimates one is usually interested to see how the solution varies for a small variation of the 
input data. Therefore, following (|5.43p . (|5.44p and (|5.45p . we assume that in the case of IPS. 2 

\\v,\\c.(s^-^<Se^^yT>0, (5.46) 

and in the case of IP5.3 

||^2llc^(P,) <'5e'*'^,Vr>0, (5.47) 

where 6 G (0, 1) is a sufficiently small number. Note that it is not necessary that S = B, where B is the number 
from (|5.43p . Indeed, while the number B is involved in the estimate of the norm ||w||(-i4(]g„^jg j^]-, , VT > in 
the entire space, the number (5 is a part of the estimate of the norm of the boundary data for both IP5.2 
and IP5.3. 

Remarks 5.1. 

1. The number 5 can be viewed as an upper estimate of the level of the error in the data ipi,ip2. Hence, 
Theorems 5.3 and 5.4 below address the question of estimating variations of the solution / of either IP5.2 
or IP5.3 via the upper estimate of the level of the error in the data. 

2. Since the kernel of the transform £ in (|5.48p decays rapidly with t — > oo, then the condition t G (0, oo) 
in (j5.44p . (|5.45p is not a serious restriction from the applied standpoint. In addition, if having the data in 
(|5.44p . (|5.45p only on a finite time interval t G (0,r) and knowing an upper estimate of a norm of the 
function / in (|5.42p . one can estimate the error in the integral (|5.48p when integrating over r G (T, oo) . 
Next, this error can be incorporated in the stability estimates of theorems of this section. 

5.3.2. Transformation to the parabolic case. Consider the following analog of the Reznickaya's 
transform (|3.38p 

Cg = g{t) = ^^l cxp(^-^^g{T)dT. (5.48) 



The transformation (|5.48p is valid for, e.g. all functions g £ C [0, oo) which satisfy \g {t)\ < Age*^"*, where Ag 
and kg are positive constants depending on g. It follows from (j5.43p that the solution u {x, t) of the problem 
(|5.4ip . (|5.42p satisfies this condition together with its derivatives up to the fourth order. Obviously 



d 








dt 






.V^^^'^i 4^. 



Hence, 

£ ig") ^ g (t) , V.g G [0, oo) such that g' (0) = 0. (5.49) 
Changing variables in (|5.48p t z,z :— T/2^/t, we obtain hnit_^o+ 5 (0 =5 (0) • Denote 

V rw. (5.50) 

It follows from (fOSl and (|09|) that 

V e C'2+aa+a/2 ^ [q, T]) , Va G (0, 1) , VT > 0. (5.51) 
By (|5.4ip . (|5.42p and (|5.5ip the function v {x,t) is the solution of the following parabolic Cauchy problem 

vt= Lv,xeW\t>0, (5.52) 

v{x,0) = f{x). (5.53) 

We refer here to the well known uniqueness result for the solution v € (72+04+0/2 ^ [0,T]) ,VT > of 
the problem (|5.52p . (|5.53p . see, e.g. the book of Ladyzhenskaya, Solonnikov and Uralceva jl04j . 

Below in Sections 5.3.2 and 5.3.3 we work only with the function v. Thus, we set everywhere below 
T := 1 for the sake of definiteness. Denote Si = dfl x (0, 1) , Pi ~ P x (0, 1) , 

Clp^ := (x, t) ^v\s,, Clp2 ^2 i^, t) ^ v |p, . (5.54) 

Then 

(^1 e (Si) ,7p2 e (p,) . (5.55) 

Let 

tpi {x, t) = d^v Isi , i/'2 (a^i = dxiv\p^ . (5.56) 

By Theorem 5.2 of Chapter IV of [104], (|5:43| and ([534]) - ([536]) there exist numbers C {n, L) ,C{P,L)>0 
depending only on listed parameters such that 

||'0i||ci+=,o/2(5^) < C {n,L) ||^i||c2+c,i+c/2(s^) , (5.57) 
||V'2||ci+°.»/2(Pi) <C:iP^I^)\\'f2\\c^+o,.i+o,/2(p^y (5.58) 

Wc now describe an elementary and well known procedure of finding the normal derivative of the function 
V either at Si (in the case of IP5.2) or at Pi (in the case of IP5.3). In fact, an analog of this procedure 
was described in Section 5.2 for the hyperbolic case, see (|5.32p . (|5.33p . In the case of IP5.2 we solve the 
initial boundary value problem for equation (|5.52p for {x, t) £ (M"\f2) x (0, 1) with the zero initial condition 
in R"\f2 (because of (|5.40p ) and the Dirichlet boundary condition v [5^= (^j. Then we uniquely find the 
normal derivative d^v \si= ipi- Similarly, in the case of IP5.3, we uniquely find the Neumann boundary 
condition dx^v \p-^ = '02 ■ Estimates (|5.57p . (|5.58p ensure the stability of this procedure. 

Therefore, problems IP5.2 and IP5.3 are replaced with a corresponding inverse problem for the parabolic 
PDE (|5.52p with the lateral Cauchy data (|5.54p . (|5.56p . These data are given at 5*1 for IP5.2 and at Pi for 
IP5.3. Uniqueness of the solution of each of these parabolic inverse problems follows from Theorem 2.2 and 
Remark 2.2. 

Using (|5:46| . (fOS]) . ([QQ)) and ([534l) - (f538l) . wc obtain 

!l<^illc2+c,i+c/2(5^) + ||'/'i||ci+=.o/2(Si) ^ {^,L,d)S. (5.59) 
Next, using (|5.47p instead of (j5.46p . we obtain 

||^2llc2+°.i+°/2(Pi) + ||V'2||ci+°.°/2(-pi) <C2{P,L,d)5, (5.60) 

where constants Ci (fi, L, d) , C2 (P, L,d) > depend only on listed parameters. It follows from (|5.59p that 
with a different constant C := C (Jl, L,d) > 

ll^illffi(Si) + IFi|Il2(Si) ^ (5.61) 
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5.3.3. Logarithmic stability estimates for IPS. 2 and IPS. 3. Theorem S. 4. Consider IPS. 2. Let 
Assumption 5.1 holds and conditions \5.40^ , {5.^6^ he valid. Also, assume that the upper hound F of the norm 
ll^/llL2(n) given, ||V/||^^jq^ < F. Then there exists a sufficiently small number Sq = So{L,fl) £ (0,1) 
and a constant Mi = Mi (L, $7) > 0, hoth dependent only on listed parameters, such that if the numher 6 in 
{5.^6^ is so small that CS G {0,do), then the following logarithmic stability estimate holds 

II f II < ^^'^ 

11/ llL2(f2) - 



In 



[csy 

Here C ~ C (fi, L,d) > is the number in i5.61\) . 

Consider now IPS. 3. Denote x = {x2, a;„) . Changing variables {x', t') = ^v^a;, dt^ with an appropriate 
constant 6 > and keeping the same notations for new variables for brevity, we obtain that without loss of 
generality we can assume that 

nc\xi + \x\^ <- xi>q\ . (5.62) 



4' 



Denote 



$ = 



{(.T,<):.TiG(0,l),a^=(a:2,a;3,...,a;„)e(-l,l)" \tG(0,l)}, (5.63) 
ai$ = $nP = |(x,t) : xi = 0,x e (-l,l)""\t e (0,1)| . (5.64) 

Recall that ()5.47p implies (|5.60p . Hence, assuming that (|5.47p holds and using (|5.63p . (|5.64p . we derive that 
there exists a constant C = C (L, $, d) > such that 

II^2|Ih1(9i$) + ||V'2|L,(a<i>) < C5. (5.65) 

Theorem S.5. Consider IPS. 3. Let Assumption 5.1 holds and {S.J^O^ , |5./^7| ) be valid. Also, assume 
that for a certain a G (0,1) the upper bound F of the norm \\f\\(j2+a(ji^ is given, i.e. \\f\\(j2+a(ji^ < F. 
Then there exists a sufficiently small number 5a = (5o {L,'^) G (0, 1) and a constant M2 = M2 (L,^) > 0, 
both dependent only on listed parameters, such that if the number S in |5.^7| ) is so small that CS G {0,Sq), 
then the following logarithmic stability estimate holds 

n n M2F 



I In 



CS' ^ 



Here C ^ C (L, $, d) > is the number from S5. 65\} . 

S.4. QRM in the hyperbolic case. Although the Quasi-Reversibility Method (QRM) was discussed 
in Section 2.5, it is worth to discuss it here again for the specific hyperbolic case. The reason is that we now 
have two types of stability estimates, which are different from the Holder stability estimate of Section 2.5: 
the Lipschitz stability estimates of Theorems 5.1, 5.2 and the logarithmic stability estimates of Theorems 
5.3, 5.4. We start from IP5.3. The case of IP 5.2 is not discussed here since it is similar. 

S.4.1. QRM for IPS. 3. To use embedding theorem, we work in this section only in 3d. The 2d case 
is similar. We now increase the required smoothness of the solution of the problem (|5.4ip . (|5.42p . To do 
this, we replace Assumption 5.1 with Assumption 5.2, where we use the norm ||'"|lc'8('^nx[o t]) instead of the 
'x[0T]) Assumption 5.1. 
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Assumption 5.2. Wc assume that integers k,p in (|5.38|) . (|5.40p . coefficients of tlie operator L and tlic 
initial condition / are such that there exists unique solution u £ C^^ (R" x [0,T]) ,VT > of the problem 
([OT|) . (jO^ satisfying 

hllc-(K"x[o,T]) <5e'^^>Vr>0, (5.66) 

where the constants B = B [L,B) > 0,d = d[L,^ >0 depend only from the coefficients of the operator L 
and an upper estimate B of the norm H/Upp^Qj < B. 

We assume in Section 5.4.1 that Assumption 5.2 holds. Hence, Theorem 5.2 of Chapter IV of the book 
[TUij . ([OS]) and imply that functions 

and there exists a number C3 (P, L) > depending only on listed parameters such that 

||V'2||c8+=.i+o/2(Pi) ^ C^{P-,L) ||<^2llcio+°.=+°/2(Pi) ■ (5-68) 
Let sets be the ones introduced in (fSJaj) . (|04l) . Then (|5^ and (fSJS]) imply that 

^2,^2 e i^^'M^l*) JI^2llH8.i(g^$) + 1^211^8.4(91*) ^ C'4(-P,4',i)||<?2llcio+°-5+c/2(p^) , (5.69) 

where the number C4 (P, $, L) — const. > depends only on listed parameters. 

Let the function v {x,t) be the one defined in Then v S C^o+^'^+^Z^ (R" x [0,r]) ,Vr > is the 

solution of the problem (|5.52p . (|5.53p . Denote 

r {x, t) — ^2 {^1 1) + — ^2 (^1 ^) + ^1 V^2 (2^, i) ; (5.70) 

w (x, = w (x, i) — r (x, t) , (5.71) 
p{x,t) = -{rt-Lr){x,t). (5.72) 

Let ($) £ . y 1^^^^ 0} . Then 

Vt-Lv^p (x, t) , (x, e w e i?o ($) , (5.73) 
t; I di<i^Q,Vxi 0. (5-74) 

To solve IP5.3 via the QRM, we minimize the following Tikhonov functional 

Ji (v) = \\vt - + 7 \HhH'S>) ^veH^i^), (5.75) 

where 7 > is the regularization parameter. The requirement v £ H'^ ($) is an over-smoothness. This 
condition is imposed to ensure that the function v € (<&): because of the embedding theorem. Indeed, 
we need the smoothness v £ C-^ ($) to apply theorems 2.2 and 3.1 of [HS]. We imposed in Assumption 4.2 
u e C^^ (M" X [0, T]) , VT > only to guarantee that v e ($) C (¥) . However, the author's numerical 
experience with QRM has consistently demonstrated that one can significantly relax the required smoothness 
in practical computation, see [HH 11001 llOlj . This is likely because one is not using an overly small grid step 
size in finite differences when minimizing functionals like the one in (|5.75p . Hence, one effectively works with 
a finite dimensional space with not too many dimensions. This means that one can rely in this case on the 
equivalence of all norms in finite dimensional spaces. Thus, most likely one can replace in real computations 
7lkllff4($) with7||z;||^2,i(<j,). 

Let (, ) and [, ] be scalar products in L2 {^) and H'^ ($) respectively. Let the function G Hq ($) be a 
minimizer of the functional (j5.75p . Then the variational principle implies that 

{dtUj — Lu, dtw — Lw) + 7 [u, w] = {p, Wt — Lw) , Vw G Hq ($) . 
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Lemma 5.1. For every function p (z L2 ($) and every 7 > there exists unique minimizer = 
Uj (p) £ iJg (<f>) of the functional |5. 75[ ). Furthermore, there exists a constant AI = M (L, (f>) such that the 
following estimate holds \\u^\\h4(^-) < Mi7~^^^ Ibllisl*) " 

Lemma 5.1 is an obvious analog of Lemma 2.5. The idea now is that if {x, t) G Hq ($) is the minimizer 
mentioned in Lemma 5.1, then the approximate solution of IP5.3 is 

f^{x)=u^{x,0) + r{x,0). (5.76) 

The question of convergence of minimizers of to the exact solution is more difficult than the existence 
question of Lemma 5.1. To address the question of convergence, we need first to introduce the exact solution 
as well as the error in the data, just as this is always done in the regularization theory [51 [T31 HH [701 1145j . 
also, see Theorem 2.7. We assume that there exists an "ideal" noiseless data G C^^ {Poo)- Following 
(|5.47p and (|5.66p . we assume that with a sufficiently small number 6 G (0, 1) the following estimate holds 

11^2 - ¥';ilc-(P,) < Se'^yr > 0. (5.77) 

Let ~ Cip2i the function /* satisfying (|5.40p is the solution of IP5.3 for the case of the noiseless data (p2, 
the function u* G C^^ (R" x [0,T]) ,Vr > satisfying tjbM^ is the solution of the Cauchy problem (|5^ . 
(|5.42p with / := /*. Following (|5.50p . denote v* = Cu*. Let functions 1^2, ■■■,p* have the same meaning 
as corresponding functions in (|5.67p - (|5.72p . except that they are generated by the noiseless data (p2- Then 
(1091) and fSJ?! imply that 

Ib-P*|li,(*) <C5($,i,d)<5, (5.78) 

where the constant C5 (P, $,i,d) > depends only on listed parameters. 

Theorem 5.6 establishes the convergence rate of the QRM. The proof of this theorem is using ()5.78p and 
is similar with the proof of Theorem 3.1 of [95]. 

Theorem 5.6. Let Assumption 5.2 and condition |5. 77\ l be valid. Suppose that the regularization 
parameter 7 = 7 ((5) := (5 G (0, 1) . Let the function G Hq ($) be the unique minimizer of the functional 
( [5. 751 ) (Lemma 5.1). Then there exists a number Y = Y{^,L,B,d) > and a sufficiently small number 
So = So {L, $) G (0, 1) such that if S € (0, i5o) , then the following logarithmic convergence rate takes place 



|A(.)-r|L,(^)<^=^==, (5.79) 



Y 

where the function f-f{s) (x) is defined in |5. 76\ l. In addition, for every uj G (0,aJo) there exists a number 
p = p{L, ^,B,d) G (0, 1/2) such that the following convergence rate takes place 



< Y5P, 



where the domain D1/2 = \ {x,t) : xi> 0,xi + \x\' 

5.4.2. QRM for Problem 5.1 of Section 5.1. In this section we assume that condition (|5.6p is in 
place. This enables us to use the Lipschitz stability estimate (|5.10p of Theorem 5.1. We consider the problem 
of the recovery of the function u {x,t) from conditions (|5.3|) . (|5.4p . We rewrite these conditions now as 

n+l 

An : = utt~ (x) Au - ^ bj (x, t) u^j - bo {x, t)u = f (x, t) , (x, t) € Q^, (5.80) 

u I g± = p{x,t) ,dnU |_5±=: q {x,t) , (5.81) 
where := Uf. Suppose that there exists a function F{x,t) such that 

FeH^ {Q^) ,F\s±^p{x,t),dnF \s±=q{x,t). 
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Denote w = u - F, G ^ f - AF. Let H'^ (Q±) = |c/ e (Q±) -.U \g±= dnU \g±^ o| . Using and 
(|5.8ip . we obtain 

Aw = G,{x,t)eQ$,w<EH^{Q^). (5.82) 

Thus, now we are concerned with finding the function w satisfying (j5.82p . QRM for this problem amounts 
to the minimization of the following Tikhonov functional 

H ^\\Aw^ G\\l^{Q^^ + 7 II^«IIh.(q±) ,weH^ (Q±) . (5.83) 

As usual, 7 > is the rcgularization parameter here. Let w-y £ Hq (Q^) be a minimizer of the functional 
(|5.83p . (, ) be the scalar product in L2 (Qt) [, ] be the scalar product in (Qt) ■ Then the variational 
principle implies that 

{Awy,Av) +j[w-y,v] = {G,Av),Wv e {Q$) . (5.84) 

Lemma 5.2 is again an obvious analog of Lemma 2.5. 

Lemma 5.2. Let the function G £ L2 (Qt) ■ Then for every 7 > there exists unique minimizer 
€ Hq (Qy) of the functional ^5.84^ . Furthermore, with a constant Gq — Cq (A, Qy) > depending only 
on listed parameters the following estimate holds 

G 

Theorem 5.7 (convergence). Let the domain Vl ~ {\x\ < R} . Assume that the coefficient c £ G^ (ifj 
of the operator A in i5.80\) satisfies conditions i5.6\) . i5. 8]) and other coefficients of the operator A are such 

that hj £G(QT),j £ [0, n+ I]. Let T > Rj ^fff^, where the constant r]^ ~ t]^ (R, d, a) £ (0, 1] was defined 

in Theorem 5.L Assume that there exists exact solution w* £ Hq {Q^) of the problem 115. 82\) with the exact 
data G* . Let w~f £ Hq (Q^) he the unique minimizer of the functional J 5. 83\) (Lemma 5.2). Then there 
exists a constant K — K {A, R, T, d,a) > such that 

11^7 - ^1U.(Q±) < K (||G- G*||^^(Q±) + V^lk1lH^(Q±)) • (5-85) 

Ln particular, if c (x) = 1, then it is sufficient to have T > R. Also, if \\G — G*\\j^^^q±^^ < S, where 6 £ (0, 1) 
is the level of the error in the data, and if ^ £ (O, 5^] , then \\wj — ^*\\h^(^q^) — 
Proof. We have 

{Aw* ,Av)+-i[w\v]^{G\Av)+-i [w* , w] , Vv e H^ (0±) . 
Subtract this equality from (j5.84p and denote w = Wj — w *,G = G-G*. Then w £ H§ (Q^) and 

{Aw, Av) + 7 [w, v] = (g, Av^ -j[w*,v],yv£ H^ {Q^) . 
Hence, Cauchy-Bunyakovsky inequality implies that 

II^^IIL(q^) + tII^IIh^(q±) ^ 

Applying Theorem 5.2 to we obtain ([^?55|) . □ 
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G 



(5.86) 



5.5. Published results. Stability estimates and convergent numerical methods for Problem 5.1, Prob- 
lem 5.2, IP5.1, IP5.2 and IP5.3 with an arbitrary time independent principal part of the operator L in (j5.37p 
were not obtained prior to the work |95| . The Lipschitz stability estimate (j5.10p is important for the control 
theory, since it is used for proofs of exact controllability theorems. For the first time, estimate (|5.10|) was 
proved in 1986 by Lop Fat Ho |121| for the equation uu — Au = with the aim of applying to the control 
theory. However, the method of multipliers of |121| cannot handle neither variable lower order terms of the 
operator L nor a variable coefficient c [x) . On the other hand, Carleman estimates are not sensitive to lower 
order terms of PDF operators and also can handle the case of a variable coefficient c {x) . 

For the first time, the Carleman estimate was applied to this problem by Klibanov and Malinsky jSlj . 
In [5T] Theorem 5.1 for equation (|5.3p with the lateral Cauchy data (|5.4p and c{x) = 1 was proved. Using 
(|5.10p . an analog of Theorem 5.7 was also proved in [ST]. Next, the result of [5T] was extended by Kazemi 
and Klibanov [72] and also by the author in [57] to a more general case of the hyperbolic inequality (|5.5p 
with c{x) = 1. Although in publications [72] jHT] [57] c = 1, it is clear from them that the key idea is in 
applying the Carleman estimate, while a specific form of the principal part of the hyperbolic operator is less 
important. This thought is reflected in the proof of Theorem 3.4.8 of the book of Isakov [Gj. Thus, Theorem 
5.1 for the variable coefficient c{x) satisfying an analog of ()5.7|) was obtained in section 2.4 of the book of 
Klibanov and Timonov [86] as well as in the paper of Clason and Klibanov [44] . The idea of [72] was used in 
the control theory by Lasiecka. Triggiani and Yao |107| , Lasiecka, Triggiani and Zhang [1081 11091 1110] and by 
Triggiani and Yao [146] . Isakov and Yamamoto j64j used that idea to prove a stronger version of Theorem 
5.1. 

To prove Lipschitz stability without condition (|5.6p . one can impose some conditions of Riemannian 
geometry on the principal part of the operator L in (|5.37p . see Bardos, Lebcau and Ranch [S], Lasiecka, 
Triggiani and Yao [TOT] [US], Lasiecka, Triggiani and Zhang [TUHl HnH 11101 , Romanov [T371[I3H] and Stefanov 
and Uhlmann [Hg. In [lOTj [lOl [103 (nHl [1371 [133 [146] Carleman estimates were used. 

Analogs of Theorem 5.7 about convergence of QRM were proved in Klibanov and Rakesh [83], Clason 
and Klibanov [IJ and Klibanov, Kuzhuget, Kabanikhin and Nechaev [91]. Numerical testing of QRM was 
performed in these references. This testing has consistently demonstrated a high degree of robustness. For 
example, accurate results were obtained in [91] with up to 50% noise in the data. 

For explicit formulas for the reconstruction of the function f (x) for TAT (IP5.1) in the case when in 
()5.37p L = A we refer to Finch, Patch and Rakesh [51], Finch, Haltmeier and Rakesh [5^, the review paper 
of Kuchment and Kunyansky [97] and Kunyansky [98] . These formulas lead to some stability estimates as 
well as to numerical methods with good performances. Another numerical method for TAT was proposed 
by Agranovsky and Kuchment [2]. 

6. Approximately Globally Convergent Numerical Method. The first step of the numerical 
method outlined in this section is the elimination of the unknown coefficient from the underlying PDF, 
which is the same step as in BK. In this section wc briefly outline the recently developed numerical method 
of Beilina and Klibanov referring for details to[4l[l4l[l5l[l6l[l7l[l8l[l9l[20l[2ll[92l[Ml HOQl [M] [1021 [l03] . 

Numerical tests arc not presented here, since they are published in these works. 

Even though the fleld of Inverse Problems is an applied one and even though MCIPs have been studied by 
many researchers since 1960-ics, the topic of reliable numerical methods for them is still in its infancy. This 
is because of enormous challenges one inevitably faces when trying to study this topic. Those challenges are 
caused by two factors combined: nonlinearity and ill-posedness of MCIPs. In the case of single measurement 
the third complicating factor is the minimal amount of available information. Conventional least squares 
functionals for MCIPs suffer from the phenomenon of multiple local minima and ravines. This leads to 
locally convergent numerical methods, which require a good first guess about the solution. However, the 
latter is impractical. 

In the above cited series of recent publications some properties of underlying PDE operators instead of 
least squares functionals were used. A very important feature of this numerical method is that it does not 
require any knowledge of neither the medium inside of the domain of interest nor of any point in a small 
neighborhood of the true solution. For the first time, the following two goals were simultaneously achieved 
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for MCIPs for a hyperbolic PDE with single measurement data: 

Goal 1. The development of such a numerical method, which would have a rigorous guarantee of 
obtaining at least one point in a small neighborhood of the exact solution without any advanced knowledge 
of that neighborhood. 

Goal 2. This numerical method should have a good performance on computationally simulated data. 
In addition, if experimental data are available, then this method should demonstrate a good performance on 
these data. 

A crucial requirement is to achieve both these goals simultaneously rather than just only one of them. 
Because of the above mentioned substantial challenges, it is natural to have the rigorous guarantee of Goal 1 
within the framework of a reasonable approximate mathematical model. Since convergence is guaranteed in 
the framework of that model, then we call our numerical method approximately globally convergent. In prin- 
ciple, any mathematical model can be called "approximate" . Therefore, the validity of our model is verified 
via a six-step procedure, see [Mj [21] . Basically this procedure includes the proof of a convergence theorem 
and computational results for both synthetic and experimental data. If convergence theorem satisfying Goal 
1 is proved and computational results are good ones (Goal 2), especially ones for experimental data, then 
that approximate mathematical model is proclaimed as a valid one. 

On the other hand, nothing works without such an approximate model. Indeed, the author is unaware 
about such numerical methods for MCIPs with single measurement data, which would: (1) simultaneously 
achieve Goals 1 and 2 and, at the same time, (2) would not rely on some reasonable approximations, which 
cannot be rigorously justified. 

6.1. Outline of the method. Let C be a convex bounded domain with the boundary dfl G . 
Let d = const. > 2. We assume that the coefficient c (x) satisfies the following conditions 

c(x)e[l,d], c(a:) = 1 for x e K3\r2,ce C" (R3) . (6.1) 

We assume a priori knowledge of the constant d, which amounts to the knowledge of the correctness set 
in the theory of Ill-Posed problems [T31 SHI IZDl I145j . Consider the Cauchy problem for the hyperbolic 
equation 

c(a;)wtt = Aw in X (0,cxd) , (6.2) 
u (x, 0) ~ 0, ut (x, 0) = (5 {x — xq) . (6.3) 

A similar technique can be developed for MCIPs for the parabolic equation c{x) Vt = Av + a{x)v, where 
either of coefficients c(x) or a (x) is unknown |14| . Equation ()6.2|) governs, e.g. propagation of acoustic and 
electromagnetic waves. In the acoustical case c{x) = 6~^(x), where b{x) is the sound speed. In the 2-D 
case of EM waves propagation, the dimensionless coefficient is c{x) = er{x), where er{x) is the spatially 
distributed dielectric constant of the medium. In the latter case the assumption c (x) = 1 for x G ]R^\J7 
in (|6.ip means that we have air outside the medium of interest $7. And the assumption c (x) > 1 refiects 
the fact that the dielectric constants of almost all materials exceed the one of the air. Equation (|6.2p was 
successfully used in [TH [HI [92] to work with experimental data, which are obviously in 3d. The latter was 
recently explained by Beilina in [22]. It was shown in Test 4 of [22] that the component of the electric field 
E (x, t) = {El, E2, E3) {x, t) , which was originally initialized, strongly dominates two other components. 

Multidimensional Coefficient Inverse Problem (MCIP). Assume that the coefficient c{x) of 
equation W.S\) satisfies condition h6.1\) and is unknown in the domain ^l. Determine the function c(x) for 
X € f2, assuming that the following function g (x,t) is known for a single source position xq ^ in h6. S\] 

u (x, t)=g (x, t) , V (x, i) e X (0, 00) . (6.4) 

The function g (x, t) models time dependent measurements of the wave field at the boundary of the 
domain of interest. The assumption of the infinite time interval in (|6.4p is not a restrictive one, because we 
work with the Laplace transform of the function u (x, t) , and the kernel of this transform decays rapidly as 
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t —7- oo. In this MCIP the data g {x,t) arc assumed to be known at the entire boundary dil. The case of 
backscattering data can be treated similarly, see Chapter 6 in the book [TJ as well as [HJ 11021 1103] . 
Consider the Laplace transform of the function u, 



w{x,s) = J u{x,t)e dt, for s > s = const. > 0. 





(6.5) 



We assume that the number s is sufficiently large, so that the integral (|6.5p converges absolutely and that 
the same is valid for the derivatives D^u,k — 0, 1,2. We call the parameter s pseudo frequency. It can be 
proven that 



Aw - s^c (x) w = -S {x - xo) , X e M^, 
lim w (x, s) — 0. 



(6.6) 
(6.7) 



Furthermore, for each value of s > the problem (|6.6I) . (|6.7p has unique solution of the form [Ml [2T| 

exp (— s |x — xqI) _ , 



w {x, s) 



47r |a; — xo| 



+ w{x,s) ,w & C^+" {R^) , lim w{x,s)^0. 

\x\^oo 



(6.8) 



Suppose that geodesic lines generated by the function c (x) are regular and c (x) is sufficiently smooth. Let 
T {x,xo) be the length of the geodesic line connecting points x and xq. Then 



\D':wix,s)\ 



2+a 



jjk / exp [-ST ix,xo)] 









} 








2 + Q 





, s — > oo, fc = 0, 1, 



(6.9) 



where M 



\2^^ = \\'\\c^+a(j2y fix,xo) is a certain function and f{x,Xo) ^ for x & Q. It is unclear how to 
effectively verify the regularity of geodesic lines for generic functions c (x) . Therefore, we assume below the 
asymptotic behavior (|6.9p without linking it to the regularity of geodesic lines. 
We have w{x,s) > [Hill]. Denote 



V (x, s) := 



Inw (x, s) 



Since the source xq ^ fl, then 



Av + s^ (Vu) = c{x),x e n. 



(6.10) 



Now we make the same step as the first step of BK. Differentiate both sides of (|6.10p with respect to s. Let 
q {x, s) = dsV (x, s) . Then 



V {x, s) = - 

V (x, s) = V {x,s) 



q {x, t) dr + V {x, s) , 
Inw (x, s) 



S2 



(6.11) 
(6.12) 



Here the truncation pseudo frequency s > s is a large number. We call V (x, s) the tail function. The tail 
function is unknown. By (|6.9I) 

>oo. (6.13) 



\V ix,s)\,^^ ^ O (-1) , |a^F(x,5)|2+„ = O (s-2) 



The number s is the main regularization parameter of our numerical method. In the computational practice 
s is chosen in numerical experiments. 
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Thus, we obtain from f|6.10p . (|6.1ip the fonowing nonUnear integral differential equation 



s s 

Ag-2s^Vg J Vq {x,t) dr + 2s Jvq{x,T)dT 



S 

2s^VqVV ~ 4sVy J Vq {x, t) dr + 2s {VVf ^0,xen. 



(6.14) 



It follows from (|6l| that 

q {x, s) = ip (x, s) , V (x, s) e dCl X [s,s] , (6.15) 

where ipi^^s) = s~'^dshi(p — 2s~^lin.(p and ip{x,s) is the Laplace transform (|6.5p of the function g{x,t) 
in (|6.4p . We have two unknown functions q and V in equation (|6.14l) . Therefore, to approximate both 
of them, we approximate the function q via "inner" iterations and the function V is approximated via 
"outer" iterations. Suppose for a moment that functions q and V are approximated in D, together with their 
derivatives D^q,D(^V, \f3\ < 2. Then the corresponding approximation for the target coefficient can be found 
via (|6.10p as 

c{x) = Av + s^{Vvf ,x'En, (6.16) 

where the function v is approximated via (j6.1ip . We have found in our numerical experiments that the 
optimal value of s to use in (I6.16P is s s. 

To solve the problem (I6.14p . (j6.15p . we assume that q (x, s) is a piecewise constant function with respect 
s. Hence, we assume that there exists a partition s = sn < sn-i < ■■■ < si < so = s, Si_i ~ Si ~ h of the 
interval [s, s] with a sufficiently small grid step size h such that 

q (x, s) = g„ {x) for s G (s„, s„-i], go = 0. 

We approximate the boundary condition (|6.15p as a piecewise constant function, qn {x) := ip^ {x) ,x G dft, 
where '0„ (x) is the average of the function ip {x, s) over the interval (s„, s„-i) . Next, a certain system of 
elliptic equations for functions q„ (x) is derived from (j6.14p using the s— dependent so-called "Carleman 
Weight Function" exp [A (s — s„_i)] , s £ (s„, s„_i) , where A >> 1, Xh > 1. Usually we use A = 50 in our 
computations. 

We solve elliptic Dirichlet boundary value problems for functions (/„ {x) sequentially, starting from qi {x) . 
An important new element of both Section 2.9 of the book [Tl] and the paper [3T] is the choice of the first 
tail function Vi^i (x) , see this section below. Having Vi_i {x) , we solve the Dirichlet boundary value problem 
for the function qi^i {x) . Next, using (|6.1ip with V := Vi,i and substituting it in (|6.10p . we find the first 
approximation ci^i (x) for our target coefficient c (x) . This is the inner iteration, i.e. when we work only 
inside of the domain Q. To perform the outer iteration in the entire space K"^, we solve the forward problem 
()6.2p . ()6.3p with c := ci,i, calculate the Laplace transform ()6.5|) at s := s and then obtain an update 
^1,2 (x) for the tail function, using (j6.12p . We repeat this procedure m times until stabilization occurs, thus 
getting functions ci.i, Vi.^, i € [1,to] . The number m is found in numerical experiments. Next, we set 
?i '■= Qi,rmCi := ci^m, V2.1 := Vi,m and repeat the same for functions (72,1, C2,i, V2.i, i S [1,"t-] ■ Similarly for 
functions g„,i, c„,i, 14,i, « e [l,m] ,ri e [l,N]. 

We now describe how do we find the first tail function Vi,i (x) . Again following the Tikhonov concept 
[3 [HI |49l [TQl 1145] , we assume that there exists unique exact solution c* (x) of our MCIP with the noiseless 
data g* {x, t) (|6.4p . We assume that the function c* {x) satisfies condition (|6.ip . Let w* (cc, s) be the solution 
of the problem (|6.6p - (|6.7p satisfying (|6.8p . Using (|6.12p . we define the exact tail V* {x, s) for s > s as 

v*{x,s) = '^^^^^^,ys>^. 

S'' 
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Assuming that the asymptotic behavior (|6.9p holds and using (|6.13p . we obtain 



V*{x,s) = ?^ + o(^] ,s->oo,x€n. (6.17) 



for a certain function p* (x) . We truncate the second term of this asymptotic behavior. Thus, our approxi- 
mate mathematical model consists of the following assumption. 

Assumption 6.1. There exists a function p* (x) € C^"*"" (jVj such that the exact tail function V* {x, s) 
has the form 

V*{x,s):=?^,ys>s. (6.18) 



In addition, 



s 

Since q* {x, s) = dgV* {x, s) for s > s, we derive from (|6.18p that 

g*(^,^) = -^. (6.20) 
s 

Set in (|6T4l) q := q*,V := F*,s = s and use ([elS]) and (lOO]) . Then we obtain the following Dhichlet 
boundary value problem for the function p* (x) 

Ap* = in 17, p* e C^+" {n) , (6.21) 
P*|90 = -5V*(x,s), (6.22) 

where ip* {x, s) is the exact function ip (x, s) , which corresponds to the function g* (x, t) . The approximate 
equation (|6.2ip is valid only within the framework of Assumption 6.1. Although this equation is linear, 
formula ()6.16p for the reconstruction of the target coefRcient c* is nonlinear. 

Assuming that the function ip {x, s) £ C^~^°' {d^l) , consider the solution p (x) of the following boundary 
value problem 

Ap^Oinil, pe C^+" (H) , (6.23) 
p\an = -s^yj{x,s). (6.24) 



We choose the first tail function as 



(x) := (6.25) 



By the Schauder theorem there exists unique solution p of the problem (|6.23p . ()6.24p . see the book of 
Ladyzhenskaya and Uralceva |105j for Schauder theorem. Furthermore, it follows from Schauder theorem as 
well as from (|6.21[) - (|6.25|) that with a number M ~ M (il) > the following estimate holds 

\VV^,i-VV\^^<M-sU{x,l)-r{xMc^+^ian)- (6-26) 

Therefore, our main approximations are (|6.18p and (j6.19p . These approximations mean the truncation 
of the term O (s^^) in (|6.17p . The goal of (|6.18l) and (j6.20p is to obtain the accuracy estimate (|6.26p for 
the first tail. We point out that these approximations are done only on the first iteration of our method. It 
follows from ()6.26p that we obtain an approximation located in a small neighborhood of the exact solution 
already on the first iteration of our method, as long as the error in the boundary data %p* (a:,s) is small. 
Convergence theorem 2.9.4 of [T3] and theorem 5.1 of [H] guarantee that all other solutions obtained in 
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the iterative process also provide good approximations, as long as the number of iterations is not too large. 
This means that we should develop a stopping criterion numerically. The latter was done in above cited 
publications about this method. 

Recall that the main point of any locally convergent numerical method is to get a good first guess about 
the solution. And this is exactly what our approximately globally convergent method delivers. Hence, we can 
apply the second stage of our two-stage numerical procedure to refine the solution [H [TH [TH [TTl UHl \19\ [20] . 
More precisely, a locally convergent Adaptive Finite Element Method (adaptivity) was applied. A good 
numerical performance of this two-stage numerical procedure was demonstrated in [H [T31 HH [HI [T51 [TH [ID] , 
including the most challenging case of experimental data [131 [HI [lO]- The adaptivity takes the solution 
obtained on the globally convergent stage as the first guess for further iterations. Suppose now that iterations 
of the first stage are stopped prior that stopping criterion is in place. Then the adaptivity still refines the 
solution quite well. For example, in Tests 2 and 3 on page 278 of [Tl] adaptivity has started from functions 
Cn,i (x) for those n which were smaller than the number N where the stopping criterion was achieved. 
Nevertheless, images of these tests were quite accurate ones. Furthermore, this two-stage numerical procedure 
led to accurate images from experimental data, see [18] and Chapter 5 of [M]. As to some other works on 
the adaptivity technique, see, e.g. Beilina [12l[T3] and Li, Xie and Zou jll6| . 

6.2. Published non-local numerical methods for MCIPs. In parallel with the above approxi- 
mately globally convergent method, another one was developed by the group of researchers from University 
of Texas at Arlington in collaboration with the author, see Klibanov, Su, Pantong, Shan and Liu [96], Pan- 
tong, Su, Shan, Klibanov and Liu |128j . Shan, Klibanov, Su, Pantong and Liu jl40| . Su, Shan, Liu and 
Klibanov |143j and Su, Klibanov, Liu, Lin, Pantong and Liu |144j . A globally accelerated numerical method 
for optical tomography with continuous wave source In this case the 2d MCIP for the elliptic equation 

Au — a (x) u = —S {x ~ xq) , x G M^, lim u (x, xq) = (6.27) 

was considered. Let F C (R^\ri) be a straight line. The MCIP consists in the determination of the unknown 
coefficient a (x) for a; g in (|6.27[) from the function (/? (x, xq) = u {x, Xq) \xGdn,xoer ■ This is the problem of 
the so-called Optical Diffusion Tomography with a direct application in optical imaging of strokes in brains. 
In this case xq is the position of the light source, u{x,xo) is the light intensity and the function a{x) is 
proportional to the absorption coefficient of light. The most important difference between this problem and 
the one discussed in Section 6.1 is that the asymptotic behavior of tails like the one in (|6.13p is not the 
case here. Thus, tails are treated quite differently in [96l I128[ Il40i 11431 1144] . In |144| accurate images from 
experimental data for a phantom medium were obtained. 

In works of Bikowski, Knudsen and Mueller |30], DeAngelo and Mueller [46], Hamilton, Herrera, Mueller 
and Von Herrmann [56] and Siltanen, Mueller and Isaacson |141| non-local reconstruction techniques are 
considered for the 2D problem of Electrical Impedance Tomography. In particular, images from experimental 
data were obtained in |46| . In papers of Alexeenko, Burov and Rumyantseva [3], Burov, Morozov and 
Rumyantseva [38] and Burov, Aleksecnko and Rumyantsevat [39] the non-local reconstruction method of 
Novikov [124l[l25l[l26l[l27] for the MCIP for an elliptic PDF with the data given in the form of the scattering 
amplitude is implemented. A non-local numerical method for a 2d MCIP for a hyperbolic PDE was developed 
by Kabanikhin and Shishlenin in [68] [69] . It is based on a 2d analog of the Gel'fand-Levitan-Krein equation. 
Just as the technique of [H] , algorithms of [3J [35] [551 [M] also use some reasonable approximate mathematical 
models, which arc not rigorously justified. 
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